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Abstract 

We construct and justify leading order weakly nonlinear geometric optics expansions for nonlinear 
hyperbolic initial value problems, including the compressible Euler equations. The technique of simul- 
taneous Picard iteration is employed to show approximate solutions tend to the exact solutions in the 
small wavelength limit. Recent work [2] by Coulombel, Gues, and Williams studied the case of reflecting 
wave trains whose expansions involve only real phases. We treat generic boundary frequencies by incor- 
porating into our expansions both real and nonreal phases. Nonreal phases introduce difficulties such as 
approximately solving complex transport equations and result in the addition of boundary layers with 
exponential decay. This also prevents us from doing an error analysis based on almost-periodic profiles 
as in [2]. 
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1 Introduction 

In this paper we considc 

(a;', Xd) — (t, y, xa) — {t, x") : a;^ > 0} for a class of equations which includes, in particular, the compressible 



In this paper we consider quasilinear hyperbolic fixed boundary problems on the domain M.^ — {x 
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Euler equations. The class consists of systems of the fohowing form: 



d 
3=0 

(1-1) / x' ■ P 

b{'"e)\xd=Q ^go + cG ( x', — — 

Ve = uq in t < 0. 

Here, Aq = I, and we assume Aj e C°°(R^,M^'), / € C°°(R^,K^), and b € C°°(K^,Ep). For the 
boundary data, we take G{x',9o) e C°°(R'^ x T^,RP) periodic in 9o and supported in {a;o > 0} and the 
boundary frequency (3 — (/3o, ■ . • , (3d^i) G M'' \ 0. 

We construct leading order weakly nonlinear geometric optics expansions and justify the expansions by 
establishing their convergence in L°° to the unique exact solutions in the large frequency limit, i.e. the small 
wavelength limit, represented by e 0. This paper completes the treatment of generic /3 in the following 
sense: we allow for all /3 but that in the glancing region, a set of measure zero which is often regarded as a 
singular case. The case of /3 belonging to the hyperbolic region, the interior of a cone determined by the Aj, 
was treated in [2] , where the authors constructed leading order expansions of highly oscillatory wavetrains 
with real phase functions. Our construction includes similar wavetrains with real phases, but in handling 
more general /3 we are required to add to the expansions exponentially decaying oscillations, which we refer 
to as elliptic boundary layers, featuring nonreal phase functions. 

1.1 The problem and its reformulations 

First we motivate and outline the relevant reformulations of the problem (jl.ip and the form of the expan- 
sion. 

1.1.1 The solution as a perturbation 

Seeking = uq + eue, a perturbation of a constant state uq with /(mq) = and b(uo) = go, we translate 
to the following problem for (where the Aj have been altered accordingly:) 

d 



(a) P{eue,dx)u^ := Aj{eue)dxjUe = T{eue)ue, 



3=0 

(1-2) / x'-P 

(6) B{eu^)ue\xd=a ^ G { x' , 



(c) = in i < 0. 

where B{v) and J-{v) denote the p x N and N x N real matrices, C°° in v, defined by 

(1.3) B{v)v = b{uo + v)-biuo), F{v)v^ f{uo + v). 
Let us define the operator 

d 

(1.4) L{dx) :=P(0,9,) = 9t + ^^j(0)9,,. 

We assume that L{dx) is hyperbolic with characteristics of constant multiplicity, making the following: 

Assumption 1.1. The matrix Aq = /. For an open neighborhood O of £ M^, there exists an integer 
q > 1, some real functions Ai, . . . , Aq that are C°° on O x \ {0}, homogeneous of degree 1, and analytic 
in ^, and there exist some positive integers vi, . . . ,Vq such that: 

d q 

(1.5) det[r/ + ^0^,(")] -n(^ + ^'^-("'^))'^' 

3=1 k=l 
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for u E O , £, = (^1, . . . ,^d) \ {0}. Moreover the eigenvalues Ai(u,f), . . . , Ag(u,^) are semi-simple (their 
algebraic multiplicity equals their geometric multiplicity) and satisfy Xi{u,^) < • • • < Xq{u,^) for all u E O, 
C€E]R^\{0}. 

Also, for our study we consider a noncharacteristic boundary {xd = 0}: 

Assumption 1.2. For u E O the matrix Ad{u) is invertible and the matrix B{u) has maximal rank, its rank 
p being equal to the number of positive eigenvalues of Ad(u) (counted with their multiplicity). 

We now provide some notation regarding frequencies r — 17 € C and rj e M.'^^^ dual to the variables t 
and y, respectively. We define the matrix 

(1.6) A{O:=-iA^\0)(^Tl + J2r^,A,{0)^ , C := (r - 17, ??) e C x M^'-\ 

and the sets of frequencies 

S := {(r - i7, 77) e C X R^-' \ (0, 0) : 7 > o} , E := {c G S : + 7^ + \r^\^ = l} , 

So := {(r, 7?) G M X M'^-i \ (0, 0)} = S n {7 = 0} , So := S n So . 

With these in mind, we define the symbol 

d 

(1.7) L(T,O:=r/ + ^C,A,(0). 

The following is a result due to Kreiss [S] in the case of strict hyperbolicity, i.e. when all the eigenvalues 
in Assumption 1 1 . 1 1 have multiplicity Vj = 1, and to Metivier [S] in our more general case. 

Proposition 1.3 ([51 |S]). Let Assumptions {T7l\ and \1.2\ be satisfied. Then for all C, E 'E,\ So, the matrix 
,4(0 has no purely imaginary eigenvalue and its stable subspace ¥f{C) has dimension p. Furthermore, W 
defines an analytic vector bundle over S \ So that can be extended as a continuous vector bundle over S. 

For (r, 77) e So, we define E*^(t, 77) to be the continuous extension obtained in Proposition 11.31 of to 

Now we describe our assumption of uniform stability, defined as in [5] [3] : 

Definition 1.4. The problem (|1.2p is uniformly stable at u = Q if the linearized operators {L(dx), B(fS)) at 
u — are such that 

(1.8) 5(0) ; E'*(r — 17, 77) — > C is an isomorphism for all (r — 17, rj) G S. 
Assumption 1.5. The problem (|1.2p is uniformly stable at u = 0. 

We now discuss an important example, the Euler equations, which satisfies the above assumptions so that 
our main results may be applied. We remark, however, that we also require Assumption 11.71 an assumption 
on the boundary frequency j3, which is explained in Section [1.1.21 

Example 1.6 (Euler equations). The following are the isentropic, compressible Euler equations in three 
space dimensions on the half space {xg, > 0}, in the unknowns density p and velocity u = {ui,U2,u^)'. 



(1.9) dt 
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where p(p) is the pressure. The hyperbolicity assumption, Assumvtion ll. 1[ is satisfied in the region of state 
space where p> 0, —p'{p) > 0. In this case the eigenvalues Xk{p,u,^) are 

(1.10) Ai - w-e-clCI, Aa ^M-e, A3 = u-e + c|e|, with (1^1,1^2,1^3) = (1,2,1). 

For this problem, the boundary condition we impose is the natural ''residual boundary condition, " which is 
obtained in the vanishing viscosity limit of the compressible Navier-Stokes equations with Dirichlet boundary 
conditions. Fixing any constant state [p, u) with U3 ^ {0, — c, c} about which to linearize the problem, so that 
{p,u) = uo in p.ip . we have noncharacteristic boundary {xj ~ 0} for the system. Consider in particular 
Assumption \1.5\ in each of the following cases: 

(a) Subsonic outflow: U3 < 0, jusj < c. In this case there is exactly one positive eigenvalue of 
A^{p,u) (p — 1), so we need one scalar boundary condition. Taking b{p,u) — in (jl.ip . we have B{0) = 
[0001] and Assumption ] 1 . 5\ is satisfied. 

(b) Subsonic inflow: < U3 < c. For this we get p — 3 and boundary condition b{p,u) — {pu^,ui,U2) 
and linearized operator 

(1.11) B{Qi){p,u) = {pu^ + pu3,ui,U2), 

which satisfies Assumvtion \1.5\ 

(c) Supersonic inflow: < c < U3. This case is trivial, with p = A and B{0) the 4x4 identity matrix, 
and so Assumvtion \1.5\ holds. 

(d) Supersonic outflow: U3 < 0, \u3\ > c. This is another trivial case, with p — 0, where B{0) is 
absent, meaning Assumption \1.5\ holds vacuously. 

With clear modifications of the above statements, the same holds for the 2D Euler equations, which are 
in fact strictly hyperbolic. For complete proofs and discussion verifying these cases, we refer the reader to 
Section 5. 

To study geometric optics for nonlinear problems, it is important to establish the existence of exact 
solutions on a fixed time interval independent of the wavelength e. For the system (|1.2p and thus (|l.ip . this 
was proved by M. Williams in |llj through the use of the singular system, which we discuss further in section 



1.1.2 The approximate solution constructed in terms of an M-periodic function 

The goal of this paper is to obtain qualitative information about the exact solution to (|1.2p by explic- 
itly constructing an approximate solution with a geometric optics expansion that exhibits the qualitative 
information, and then showing that this approximate solution tends to the exact solution as e — > 0. 

Regular boundary frequencies. 

We now discuss an assumption regarding the boundary frequency /? = (t, r/) e M x K'^^^ \ (0,0) which 
is relevant to the construction of the approximate solution. The planar real phase (po is defined on the 
boundary by 

(1.12) Mx')-=Tt + ri-y^(3-x'. 

Note the relationship to the boundary data (jl.2p fb). The oscillations on the boundary associated to 4>o result 
in osciUations associated to some planar phases which are characteristic for L(d^) and which have trace 
on the boundary equal to (pQ. Such oscillations largely motivate the ansatz for our approximate solution, 
dependent on the planar phases The following assumption ahows us to obtain the 

Assumption 1.7. For = {t ~ i'j.rj) E C x M."^^^ consider the matrix 

(1.13) lA{T-ij,rj) = -A^\0) |^(r-i7)/ + E^,-4,(0)j . 

Let w^, TO — 1, . . . , M denote the distinct eigenvalues of i-4(r, rj). We suppose that each uj_^ is a semisimple 
eigenvalue with multiplicity denoted by fim ■ Moreover, we assume there is a conic neighborhood O of j3 in 
C X \ {0} on which the eigenvalues of i^(C) «re semisimple and given by smooth functions uimiOj 

TO = 1, . . . , M where — LOmiP) and uj„i{C,) is of constant multiplicity p„i- 



4 



Remark 1.8. Using Assumption 1 1.7[ one can show that the smooth functions oJmiO in fact analytic on 
a conic neighborhood of /3. In particular, the analyticity is used to prove Lemma 1 1.1 31 

We call /? a regular boundary frequency provided Assumption 11.71 holds. Now, we define the phase 
functions 

(1.14) 0m (a;) := M^') +^„,a;rf = ■ x, 

and we have the complex characteristic vector field associated to (jfm- 

(1-15) := d,,+J2-d^,^mmd,,, 

j=o 

which is real for real a;„. Moreover, for each which is real. Assumption 11.11 implies there is a unique 
km S {1: ■ • ■ ; such that r + Afc^ (?], oj^) = 0, and in this case we have — With this in mind, we 
make the following definition: 

Definition 1.9. (i) Form such that uj,-^ is real, we call (/3,w„j) a hyperbolic mode if 

(1-16) %A,„(^,w,J^O. 

(ii) For m with nonrealu^, we call an elliptic mode. 

Remark 1.10. (i) For each real w^. Assumption 11.71 guarantees that (/3,w„) is a hyperbolic mode. That 
(|1.16p holds is a consequence of the semisimplicity of w^H Also, the condition (| 1.1 61) and the implicit 
function theorem imply that a;„i is real and of multiplicity i>k,„ in a neighborhood of /3. Thus, there is a 
slight redundancy in Assumption 11.71 when w„j is real. 

(ii) When is nonreal, the fact that the Aj are real implies that iA{f3) also has the complex conjugate 
w„ as an eigenvalue, and thus — a;^/ for some m' ^ m. Furthermore, if the vector r S is an 
eigenvector of iA{l3) associated to w^, then r is an eigenvector associated to w„ — a;„,. 

(Hi) The boundary frequency /3 lies in the hyperbolic region (as defined in [5]) if and only if Assumption 
11.71 holds with all real. The elliptic region consists of all /? such that Assumption 11.71 holds with all 
nonreal. 

Example 1.11. We return to the 3D Euler equations, linearizing about (p, u) — (p, ui,U2, u^), as in Example 

EH 

For jual > c the set of regular boundary frequencies is all of MP \ 0. In the case [u^l < c, one finds that 
the set of regular boundary frequencies is 



(1.17) 




We remark that in the former case, the hyperbolic region is also Ti. = \ Q, and in the latter case 
H ^ |(r,?7) e : \t + uiiji +M2'72| > " "3 M]- 

For real we have the associated real group velocity: 

(1-18) v,„ :- VAfe„(?7,^„J, 

which has the following relationship with the characteristic vector field (jl.lSp : 

1 df Xh (ri,uj^) 

Observe that each group velocity v™ can be thought of as either incoming or outgoing with respect to the 
interior of the domain W!^, in particular since the last coordinate of v.,„ is nonzero, by (|1.16p . With this in 
mind, we classify the phases in the following way: 

^ Assumption [LTI has no immediate implication for the multiphcity of nonreal 

^If u]_^ is real and d^^Xkm iVt^m) ~ ^' refer to {P,Ld^) as a glancing mode. An explanation of how III.I6I I follows from 
semisimplicity can be found in the proof of Lemma 2.7 of [§]. 
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Definition 1.12. For real lo_^, the phase (j)m is incoming if the group velocity is incoming (that is, 
d^d^k„S'nT'^m) > Oj; it is outgoing if the group velocity v,„ is outgoing ('S^dAfe^ (?7,a;,„) < 0). 

With this classification of the real phases, we let I denote the set of indices m S {1, . . . , M} such that 
(pm is incoming and O the set of m such that 0^ is outgoing. We classify the remaining complex phases (pm , 
which correspond to nonreal a;,„, by the distinction that V is the set of m such that Im a;„j > and Af is 
the set of m such that Im w„j < 0. We thus form the partition 

(1.20) {1,...,M} ^lUOUT'UTV. 

Recall the uj,n{T,r]) are eigenvalues of i^(r, 77), meaning the iuj„i{T,ri) are the eigenvalues of A{T,ri). 
Consider also the fact that Im uim{T,ri) > if and only if Re iuJm{T,ri) < 0. From these observations, we 
see that the stable subspace (r, ij) of A{t, rj) must contain each of the eigenspaces corresponding to some 
uJm{T,ri) with Re «a;„ < 0, so that E*(r, 77) contains the subspaces Ker L{d(j)m) for m e V. The incoming 
phases, corresponding to m g I, also play a role in setting up a decomposition for the stable subspace E" at 
the boundary, which is established in the following lemma. 

Lemma 1.13. The stable subspace E''(r, 77) admits the decomposition: 

(1.21) E^(t,^) = e^eiupKer L(d0™), 

where, in the decomposition, the vector spaces Ker L{d(j)m) for to e I admit a basis of real vectors. 

Proof. It is easy to show that the subspaces Ker L{d(j)rn) for m G V are in E^(t, 77). We will show that this 
is also the case for the subspaces Ker L{d(j)m) with m Gl, from which the result will follow. Since E^(r, 77) 
is close to E'*(t — 17, 77) for small 7 > 0, in accordance with Proposition 1 1.31 it will suffice to show for to G I 
that ii^mix. ~ ili'n)i close to iw„j, satisfies 

(1.22) Re za;m(T- 77,77) < 0. 
Recalling (I1.19p . we see that since (pm is incoming, i.e. to G I, we have 

(1.23) 5cot^™(r,7/) < 0. 

Using analyticity of uJmiO with ()1.23|) . since uJmizTV) h^-s imaginary part equal to zero, it follows that 
ujjri{z~ iliV) h^'S positive imaginary part for small positive 7. Therefore (jl.22p holds. The statement that 
each of the vector spaces Ker L{d(l>rn) for to G I admits a basis of real vectors follows from Lemma [1.14[ 
which was proved in [1]. □ 

The ansatz. 

The following lemma, proved in [T], provides a decomposition which is key to the construction of our 
ansatz and will serve us later in the construction of the projectors which give us the profile equations. 

Lemma 1.14. The space admits the decomposition: 

(1.24) C~ = ®*f^iKeri(d(/),„) 

and each vector .space in (|1.24p with m G TUO admits a basis of real vectors. If we let Pi, ... , Pm denote the 
projectors associated with the decomposition (jl.24l) . then for all m — 1, . . . , M , there holds Im A^^L{d(f>m) = 
KerP™. 

For each to = 1, ... , M, we choose a basis of Ker L{d(pm). 

(1.25) rm,k, fc = 1,...,^. 



m 1 



where for real uj_^, we take a basis of real rm,k- Now we are in a position to give the ansatz for our approximate 
solution to (|1.2p . 

(1.26) u'iix) = v{x) f-^' ^ 

m=l fe=l ^ 
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where v{x) and the (Jm,k{x,9m) are functions. Additionally, we require that each of the (Jm,kix,Om) is 
periodic in 0m with mean 0, and we will refer to these as the periodic profiles. With these, we plug (|1.26p 
into (|1.2p (a) and expand the error in ascending powers of e: 

(1.27) P(e<, - J-(eO< - e"' ( J] Xi I^id<Pm)de^cT,n,k U, J + e" (...) + ... . 

Observe that since each r,n,k belongs to Ker L{d4)„i), the term of order i in (|1.27p is in fact zero. 
Equation (|1.26p is a special case of a substitution of the form 

(1.28) <(.)=VO('x,^^(") 



e e 

for a function V°(x, 9i, . . . , 9m) periodic in . . . , 9m)- One finds that solving (ll.2p (a') to order - amounts 
to satisfying the condition 

M 

(1.29) C{dg)V'' := J2 Hdcj)m)dg„y = 0, 

m— 1 

where L{d(j)m) — A'^^ L{d(j}m) ■ Indeed, (|1.29p is satisfied if V° has the form 

M Mm 

(1.30) V°ix,9i,...,9M) = v{x) + J2Y. 

^ra,k {x ^ 9m)^m,k : 

m=l k=l 

which is consistent with the ansatz of (|1.26p . 

One can try to solve (|1.2p (a) to higher order by constructing a corrected approximate solution which, 
upon replacing in ()1.27p . results in the vanishing of the coefficients of higher powers of e. In fact, we will 
obtain more conditions on v and the <7m,k which will ultimately determine our ansatz for by considering 
the existence of such a which corrects and agrees in the leading term. Agreement is meant in the sense 
that, for some Af-periodic V^{x,9i, . . . , 9m), we have 

(1.31) 

(f>iix) 4'Mix)\ g/ (j)i{x) (I)m{x)\ , ..if (t)i{x) (t>M{x) 



<(x) = <(x)+.Vi x,^,...,^^^^ =V° x,^,...,^^^^ 



x, 



Here is not generally expected to share the form of V° seen in (|1.30p . as we do not expect to solve 
()1.29p . The function is determined by other differential equations arising from setting higher order terms 
of (fr27l) to zero. 

Remark 1.15. Ensuring the possibility of constructing such w^, which suggests we have the right leading 
term requires solvability conditions to hold; these will be translated into conditions known as the profile 
equations on the profiles, that is, vix) and the crm,k{x,9m)- We will discuss the profile equations in greater 
depth later. 

Expansions with elliptic boundary layers were constructed in |10) . treating the semilinear problem for 
generic and in [5], where a semilinear dispersive problem with maximally dissipative boundary conditions 
was considered. In both papers, higher order expansions in ascending powers of e, such as 

(1.32) V \ x, +... + e V \x, 



are used to justify the approximate solution as well as construct the exact solution. High order expansions 
involving surface waves, which arise from a failure of the uniform Lopatinski condition at a frequency /3 
in the elliptic region, were constructed and justified in 7 for quasilinear hyperbolic systems. Here in the 
spirit of [5] we justify a leading term expansion for solutions to (jl.ll) oscillating with multiple phases. In 
this situation it is impossible to construct a high order expansion without a small divisor assumption, an 
assumption we do not make. While the small divisor assumption would exclude /? only from a set of measure 
zero, verifying it for a given /? can be difficult if not impossible. 

^In fact, lO^ also treats /3 resulting in glancing modes of order two, which do not result in blow-up, unlike higher-order 
glancing cases. 
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1.1.3 The solution of the singular system and the main theorem 

The following theorem verifies that the approximate solution is in fact close to the exact solution for small 

e. 

Theorem 1.16. Suppose the assumptions hold that L{dx) is hyperbolic with characteristics of constant mul- 
tiplicity, the boundary {xd = 0} is noncharacteristic, {L{dx), B{0)) satisfies the uniform stability condition, 
and P is a regular boundary frequency (i.e. Assumptions HOI \1.2l \1.5l resp.) Then given the exact 

solution Me of ( |1.2[) defined on a time interval {—oo,T] independent of e, for as in (jl.26p . where v and 
the (Jra,k satisf^ the profile equations, defined on a time interval (— oo,T'] independent of e, we have 

(1.33) \u,{x) - <(a^)Uoo((_^^^„],s:t) ^ «^ ^ ^ 0' 

where Tq is the minimum of T and T' . 

Observe for this result that we need in particular a positive lower bound on the existence time of the 
exact solution Ut to (|1.2|) which holds uniformly for small e. That is, for some positive eo, 



(1.34) inf T, > 0. 

(0,eo] 

However, applying the standard theory to the problem (|1.2|) yields existence times Te which shrink to zero 
as the Sobolev norms of the initial data blow up in the limit e — )■ 0. To get estimates uniform in e, one may 
use a reformulation of the system (|1.2[) known as the singular system, which is also used in the framework 
of [5]. In fact, to prove the main theorem, we show a stronger result involving the solution of the singular 
system. The singular system is obtained by recasting (11.21) in terms of an unknown J7e(x, 6q) periodic in 9q 
such that a solution of (|1.2p is to be formed by making the substitution 



B ■ x' 

(1.35) u^{x)^UAx,-— 



The singular system is written in the form 

(a) dx,U, + i,(eC/,) (dx^ + ^ j - F{eU,)U,, 
(1.36) 

(b) BieU,)iU,)U,=o^G{x',9o), 

(c) Ue=Omt<Q, 

where F = A^^T. Indeed, the main benefit of expressing the problem in terms of the new unknown Ue 
is that it is possible to obtain estimates uniform in e for Ue and thus existence of the solution Uc to (|1.2p 
in a fixed time interval independent of e, demonstrated in |11| by M. Williams, where he constructed the 
exact solution and proved these estimates which are crucial to our analysis as well as the analysis of [2]. 
The solution Uf: of the singular system is constructed in Theorem 7.1 of [TT], with the use of the following 
iteration scheme: 

(1.37) 

&)i?(ef/r)(C/f+i)U,.o = G(x',0o), 
c) U'^+^ =Omt<0. 

Once the are obtained, C/e is found by taking the limit as n — cxd of the U^. 
The stronger result we will show in order to get Theorem II. 161 is the following: 



*The profile equations are not solved exactly. Description of the sense in which error terms for the profile equations are to 
be small can be found in the error analysis. Essentially, the error must satisfy the hypotheses of Proposition 12.40] 
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Theorem 1.17. Define 

(1.38) U^ix, eo) = V" (x, 00 + u^i^, . . . , 00 + ^My) 

where V° is as in (jl.30l) . If G{x' , 9q) e for sufficiently large s, then for the exact solution of p.35p . 
we have 

(1.39) |C/,(a;,0o)-Z^°(a;,0o)|L~(x,.ff»-i(.',e„)) ^0 ase^O. 



To handle the supremum norms taken over Xd € [0, cjo) and to accommodate the appearances of ^ in 
the M periodic arguments of (|1.38p . we introduce the placeholder = Thus, we consider a function 

(1.40) U\x, 00, £.d) ■■= V°(x, 6*0 + WiCrf, . . . , 00 + (HmU) = v{x) + '^"^A^^ + Olm£.d)rm,k, 

ni—1 k—1 

which we may also write as 

M Mm 

(1-41) U°{x,9o,U) + Xm^":'^(2;,0o,Cd)r-™,fe, 

m— 1 k—1 

where profiles 4'm,kix,9o,^d) are given in terms of the periodic profiles (Jm.k{x,9m) by 
(1-42) ■0m,fe(a;,0o,C<i) := (T„,fe(a;, 6*0 + ^iimCd)- 

As a result, assertions such as Lemma [2.421 and Proposition 12.401 which concern substitutions of the form 
(1.43) U^A^.Oo) =Z^°(x,0o,e<i)lc,=^, 

play an important role in the error analysis. 

Remark 1.18. When a;,„ is real, since (Tm,k{x,9m) is periodic in 9m, the profile 'ipm,kix,0o,^d) defined in 
(|1.42p is almost-periodic in (00, ^d)- Furthermore, if /3 is in the hyperbolic region defined in Remark ll.lOf iii). 
all the are real, and so in that case U'^{x,9o,^d) is almost-periodic in (0o,^d). In the case that some of 
the are nonreal, in order to make sense of the substitutions (ll.42p we must first extend the (Jm,k{x, 9m) 
to be defined for 0,„ in the complex plane. This process is detailed in Section ri.2.11 

1.2 Role of nonreal phases and the resulting boundary layers 

We now point out a key feature distinguishing our study from that in ^ . In it is assumed that /3 belongs 
to the hyperbolic region. This is equivalent to requiring that both (i) /? is a regular boundary frequency (i.e. 
/? satisfies Assumption II. 7p and (ii) all the eigenvalues w„ of the matrix jA{f3) are real. In our study, we 
allow for any regular boundary frequency /3 and must thus handle cases in which some of the cj„ are nonreal. 
Hence, functions of complex variables must be considered in, for example, (|1.26p . p.BOp . ()1.38p . and ()1.40p . 

1.2.1 Hyperbolic and elliptic profiles, and the emergence of the elliptic boundary layer 

For each m such that w„j is real, we call the <Jm,k{x, 9m), k = 1, . . . , fJ-m, hyperbolic profiles, and we note 
that in order to make the substitution 0„j = 0o + uim^d, such as those made in ()1.40p . a hyperbolic profile 
needs only to be defined for real 0,„. On the other hand, when uJm is nonreal, we call the <Jm.k{x, 9m) elliptic 
profiles. For each elliptic profile, as we vary the parameters 0o G M, > 0, the value 9m = 0o +^m^<i varies 
throughout one of the half complex planes {Im 9m > 0}, {Im 9m < 0}, depending on the sign of Im a;„. To 
make sense of substituting 9m = 0o ^m^d i'^to the argument of an elliptic profile, we define the profile first 
for real 9m and then holomorphically extend the 0„i-dependence into the appropriate half complex plane. 
To do so, we make use of the Fourier expansions of the profiles. 
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Consider a profile <Jm,k{x, 0„i), periodic in 6m with mean 0, and its expansion of the form 

(1.44) a^Ax,Om)^ J2 ajix)e'''-. 

For the moment we proceed formally, assuming that one can evaluate the above expression at 6„i = Oo+ui^S^d 
by substituting 9m = Oo + ^m^d in the argument of each exponential in the expansion for (Jm.k{x, 0m) and 
yield a convergent expansion for the profile "ipm.kix, Oq, ^d) as defined in (|1.42l) . That is, given (|1.44p . we also 
have 

(1.45) ^mA^,eo,^d)^ J2 a,(x)e^^'(«"+^™«^). 

We remark that, certainly, we have yet to make sense of this sum in the case that e C \ K, but also that 
the space of convergence must be clarified when uj^ € M. The convergence is made rigorous for both cases 
with Proposition 12.391 Let us rewrite such an expansion in the following form: 

(1.46) iJmA^. ^0, (d) = fy^nA^, Oo + uimU) = E ix)e'^'>" 6^^''^^^^^^^' c-'"^^^^^^^' . 

jez\o 

Suppose (Tm,k{x, 9m) is elliptic. We first consider the case with Im Lo_m > 0- Then, the terms in the sum in 
()1.46p with J < grow exponentially with £^d- It is easy to check that such terms result in an unsatisfactory 
candidate for our approximate solution u'^{x), as defined in (|1.26p . which is nonphysical in the sense that it 
blows up in L°° as e — > 0. This leads us to construct <JmA^^ ^m) such that aj{x) = for j < 0, and so an 
elliptic profile with Im > is to have in fact the following expansion rather than the one in (jl.44p : 

(1.47) (TmA^^^m)^ aj(2;)e*J"™, for Im > 0. 

Provided the above is a Fourier series in real 9m converging in H^'^^(x, 9m), one can show that it holomor- 
phically extends into {Im 9m > 0}. Now, for the corresponding profile ipm,kix,9o,£^d) = o'm,fc(a;,6'o +w„^t;), 
we take (|1.45p and (jl.46l) . apply the fact that aj{x) = for j < 0, and get 

(1.48) 

V'm,fc(a;,6'o,e<i) - J2 aj(a;)e*^'(^°+^"«'') = ^ aJ-(x)e*J■^°e'^'«■'=(^^-)^''e~^■^"^(^^'")^^ for Im Ci;„^ > 0. 
jez+\o iez+\o 

As a result, the profile i'm,k{x,9o,^d) must decay exponentially in £^d- 

In the case that Im lo^ < 0, similar considerations lead us to construct a corresponding elliptic profile of 
the form 

(1.49) crm,fe(2:,6'„) = aj(a;)e^J"", for Im a;„ < 0, 

a sum which holomorphically extends into {Im 9m < 0}. This kind of elliptic profile also results in 
'0m,fc(a;, 6*0, '^d) — crm,kix,9a +^j^^rf) which decays exponentially as increases. 

For a hyperbolic profile crm,k{x,9m), which has Im w„ — 0, we do not make such restrictions on its 
coefficients aj{x), so we describe it with (|1.44p and ipm.k{x,9o,^d) with the expansions (|1.45p and (|1.46p . 

With the following remark, we summarize the different forms of expansions for profiles. 

Remark 1.19. For Zm = {n e Z : nlmw,^ > 0}, i.e. 

r Z if a;,„ G M, 

(1.50) Zm = < Z+ if Im cj„^ > 0, 

[Z^ if Im < 0, 
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each periodic profile has an expansion of the form 

(1.51) cr„,fe(a:,0m) = aj(a;)e'-'^", 

and each profile ipm,,kix,9o,^d) has 

(1.52) i^mM^,eo,^d)^ J2 a,ix)e''^'°+^^^'^ ^ ^ a,(x)e^^'^«e^^'^'=(^")«^e-^'i'"(^^")«''. 

Since the elliptic profiles (ym.k{x,Om) result in ipm,kix,OQ,£,d) which decay in f^, upon replacing the 
placeholder with ^ we see they contribute to a boundary layer with rapid exponential decay in Xd for 
small e, i.e. the elliptic boundary layer, appearing in the approximate solution described by (ll.26p . 

It appears that this paper is the first to rigorously justify leading order expansions involving multiple real 
and nonreal phase functions (f>m{x), and thus both hyperbolic and elliptic profiles, in quasilinear hyperbolic 
boundary problems. 



1.2.2 Loss of almost-periodicity and construction of the corrector 

In [2], the authors constructed leading order expansions similar to ours, though aspects of the analysis 
do not work with nonreal phases. As discussed in Remark 11.181 in the case studied in [2 , all the are 
real, a condition which generally yields U'^{x, 9o,£,d) which is almost-periodic in (0o, £.d)- Thus, in that study, 
rather than having to define a kind of convergence of infinite sums of the form (ll.46p . since the functions of 
interest were all almost-periodic, it sufficed to work in the space V^^, defined by taking the closure of the set 
of finite trigonometric polynomials in defined by 

(1.53) £^^{Uix,eo,U) ■■ sup\Ui-,-,U)\E^ <oo}, 

id>o 

where 

(1.54) = C{xd, H^ix\ 0o)) n L^ixd,H^+\x', 0o)). 

However, in the case that some of the are nonreal, one loses the almost-periodicity of U^{x, 9o,£,d)- The 
introduction of exponential decay in £^d of the corresponding ipm,kix,0o,id) presents serious obstacles to 
straightforwardly adopting the approach of |2] . Many estimates proved in [2] appear to have no analogue in 
our study. 

As an example, we consider the construction of a projector used in [2], which we denote by P. The 
operator P is first defined on finite polynomials of the form 

(1.55) H{x,0o,U)^ E ffK(x)e*««^«+*'^''«^ 

and is constructed such that (PiJ)(x, ^o, ^d) = if and only if there exists a solution U of 

(1.56) C\deo,d^,)U = H{x,0o,Cd), 

where C'{dea,d^^) = A{(3)dgg +d^^. Now suppose that, in an attempt to account for the existence of nonreal 
eigenvalues aj„j of jA{l3), we take iJ to be a finite polynomial of the form 

(1.57) H{x,0o,id)^ Yl i^K(x)e^«°''°+"'-^^ 

K;=(K;o,Kd)6Zx{zGC:Im 2>0} 

noting the decay of H in ^d- Given the existence of more than two nonreal eigenvalues of jA{f3), we found 
that a linear operator P satisfying the same for such H in any sensible function spac^ typically is unbounded. 

^The nonlinearity of HI. Ill forces us to consider the action of C'(dgg,d^^) (and thus P) on vectors such as re*'^<i^<i € for 
all Kjj in some set closed under Z-linear combinations; in the event that more than two elliptic modes exist which are pairwise 
independent over Q, it follows there exists a sequence re"'d^d in the kernel of P converging in L°° to an element not in the 
kernel. 
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On the other hand, the assumptions in [2] lead to having no such eigenvalues, and in that case the projector 
P can indeed be continuously extended to P^. In fact, F : — ?> is used in [2] to obtain a corrector term 
for the expansion which is crucial to the error analysis of [Sj- We also need a corrector but we construct it 
in a different way. 

To obtain a satisfactory corrector term, we found there was no need to define a projector on the space 
of almost-periodic profiles, such as P. In fact, we were able to dispense entirely with almost-periodic for- 
mulation of the profile equations, such as the one used in [2]. Instead, it was sufficient to solve profile 
equations and use projectors in just the space for M -periodic functions, which contains the V*'(x, 0i, . . . , 0m) 

and V^a:, 6*1, . . . , 9m) featured in and f^T^ . namely H^iW^ x T^^, and therein we found we could 

naturally phrase much of the error analysis and construct a useful corrector term. In much of the error 
analysis, in place of the almost-periodic profiles with expansions in terms similar to (|1.55p . we deal with cor- 
responding Fourier expansions in terms e*"'^, a e (see Definition l2.1[ ) 6 — {6i, . . . , 9m), holomorphically 
extended in to a selected subset C^^ of C^^ (see Remark l2.3[ ) circumventing the use of almost-periodicity 
of [2]. The key observation was that for the construction of an appropriate correcter in iJ'*(M'|_^^ x T*^), 
hyperbolic objects as well as elliptic objects could be handled simultaneously, as the terms e*"'^ behave no 
worse for 6* € C*^ C than they do at e M*^. 

Interestingly, working with projectors only in i/*(M^ x T*^) to obtain the corrector appeared to be 
significantly simpler than attempting to do so in a space similar to the one used in [5]. For example, 
the continuity of our projectors is almost immediate |j while even in the case treated in |2], the proof of 
continuity of the projector P takes some effort. The trade-off for the gained simplicity is that we must 
describe convergence of expansions such as (|1.46p carefully with Proposition 12.391 instead of just using a 
space which is a closure of trigonometric polynomials and assuming all the functions needed are almost- 
periodic. Finally, while we do as much of the error analysis in ^^''(IR!^^^ x T*^) as we can, we are forced to 
introduce the substitutions 9 — {9q + c^ii^d, ■ ■ ■ ,Oo + (ilM^d) and £,d — ^ use [IT for the error analysis, 
with Proposition 12.431 

1.2.3 Complex transport equations and resonances 

In contrast with the real transport equations of the profile equations for the hyperbolic profiles, the complex 
transport equations for the elliptic profiles are generally not solvable. In |10j . a Taylor expansion in Xd is 
developed to approximately solve the corresponding complex transport equations of the profile equations, 
but using the usual Taylor error bounds would require us to consider C'^ regularity in place of regularity, 
presenting difficulties to our strategy for leading-order justification. Interestingly, however, in our proof of 
Theorem 12.411 we were able to show that merely requiring the complex transport equations of the profile 
equations to hold to first order at Xd — is sufhcient for achieving L°° convergence. The error from the 
complex transport equations belongs to a class of functions in Hj^(]S.^ , T^^) which are elliptically polarized^ 
This error term is thus handled by an application of Proposition l2.40l which clarifies the sense in which such 
functions are small if they are zero at the boundary {xd — 0}. 

Another unusual feature involving elliptic phases in the analysis is the handling of complex resonances, 
which result in couplings amongst the transport equations for the profiles. Careful examination of these 
resonances revealed that we could first work strictly within a subsystem of real transport equations not 
dependent on the elliptic profiles and solve for the hyperbolic profiles. On the other hand, resonances 
between hyperbolic and elliptic phases do affect the complex transport equations for the elliptic profiles. 

While, as in [2], hyperbolic profiles cr,n^k are obtained as limits of f. via an iteration scheme, there is no 
need to develop such a scheme to obtain the elliptic profiles. We construct satisfactory elliptic am,k by using 
finite propagation and regularity properties resulting from solving wave equations where time is thought of 
as the Xd variable rather than the t variable and using a trick of taking a frame of reference moving at the 
speed of propagation to ensure the am.k are supported in t > 0. The conditions which the elliptic profiles 
must satisfy at a;^ = are straightforwardly read off from linear relations with the hyperbolic profiles at 
Xd = and G. Inserting the values of the profiles at = into the complex transport equation evaluated on 
the boundary and requiring this to hold determines the values of the Xd-derivatives of the elliptic profiles on 

^See E and appearing in Definition 12.61 and Remark l2.7l 

^These functions satisfy EeV = V, where Ee is tfie elliptic projector defined in 112.2311 . 
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the boundarjo- However, we still develop a sequence of elliptic 'iterates' cr^ j.; we emphasize that while they 
do converge to the elliptic am.ki these elliptic iterates are not used in the construction of the elliptic am.k- 
The elliptic iterates are useful because they are convenient for the error analysis, allowing us to perform 
much of the analysis on the elliptic and hyperbolic parts simultaneously and in as uniform a way as possible. 

1.3 Use of simultaneous Picard iteration in the error analysis 

In the proof of Theorem 7.1 of [TT], for some Tq > 0, the iteration scheme for the singular system, (|1.37p . 
is used to produce U^{x,9q) and iterates f7"(x, ^o) bounded in E^^ uniformly with respect to n and e and 
which satisfy 

(1.58) lim [/" = in uniformly with respect to e € (0, ep], 

where solves the singular system (|1.36l) . 

In Section [^751 Proposition 12.301 of Section and Definition 12.311 we construct a function V^{x,9) E 

H'^{M.'^^ X T*^) which approximately satisfies the profile equations, and iterates V'''"(x, 0) approximately 
satisfying the equations of a corresponding iteration scheme. The iterates V°'" are bounded in H^^^ uniformly 
with respect to n and satisfy 

(1.59) lim V°^" = V° in . 

Again, we make the substitution seen in (|1.40p . plugging in 9 = [do + w^^d, . . . , + ^iZi/^d), to get 

(1.60) U^'^'ix, 9o,Cd) V°'"(a;, 9o + • • • , + ^M^d), 
followed — yielding 

(1.61) U^ '^ix, 9o) := V"-" (x, 00 + ^1^, . . . , ^0 + ) • 

Thus, by using (|1.59l) and applying the estimates given by Proposition 12 . 391 and Lemma [2.421 we get that 

(1.62) lim = U° in E^T^ uniformly with respect to e G (0, eo], 

where 14^ is as in Theorem 1 1.1 71 Therefore, to conclude lime_j.o^"(2;, 9o) — Ue{x, 6*0) = in Efj7^ and finish 
the proof of Theorem 1 1.1 71 it is sufficient to show 

(1.63) lim \U^^" - C/rifij-i = for aU n. 

Indeed, p.63p is proved in Section [2.61 bv induction on n. One might try to prove the statement in this 
way by applying the estimate for the linearized singular system of Proposition l2.43l to (U^'"^^ — U^~^^) . The 
problem with this is that if we take A{eU"), which we define to be the operator appearing in the left hand 
side of the equation for i.e. (|1.37l) (a). and apply it to the difference — U^~^^), the resulting 

quantity does not necessarily converge to zero as e — ^ 0. To illustrate this point, consider the relation of 
to V°. Recall that is constructeqfl to satisfy (|1.29p so that in (|1.27|) the order - terms are annihilated, 
noting that the terms of orders e, e^, ... shrink to zero with e, but observe that we have yet to handle the 
0(1) terms. Similarly, plugging in W^'""*"^, an approximation for U^, in place of U^~^^ in the equation for 
jjn+i j-gg^iits in an 0(1) error. 

This could be handled if a corrector for y*^'"^^ could be constructed, say V"'^, analogous to the 
appearing in (|1.3ip . which resulted in the elimination of the 0(1) terms upon the replacement of V^'""*"^ by 
yO,n+i ^ ^ similar equation. Then one obtains a corrector for U^'^~^^, the function U^{x,9o), from 

V^(x, 9) by plugging in 9 = {9q + Wifd, ■ ■ ■ ,9o + w^/Cd) and = ^ as before. In fact, applying A(eZ/^°'") to 
(-^o,«+i _^ ^1 _ jjn+i^ results in an error which converges to zero as e — ^ given the existence of suitable 

®For further discussion, see Section 12.41 

^Strictly speaking, details on the construction of V'' come later in the paper with the solution of the profile equations. 
Instead, this refers to imposing the condition that V" has the form l|1.30|l . 
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V"'^. After an application of the estimate from [TT], (|1.63p can then be obtained, assuming dA} is bounded 
by Ce in an appropriate norm. 

However, there are two major obstacles to constructing such a corrector for V°'"+^. The first is 
that each of our iterates V"'""*"^ only approximately solves the corresponding iterate equation of the profile 
equations. This is forced on us since our profile equations include complex transport equations which are not 
generally exactly solvable. On the other hand, as mentioned in Remark II. 15[ typically in geometric optics, 
the solvability conditions on an ansatz such as V° which imply V° has a corrector are satisfied by exactly 
solving profile equations; similarly, in [5], solvability is achieved through exact solutions of the iterated profile 
equations. Instead, however, we found it was sufficient to solve each of the iterated profile equations up to 
an error which is zero at the boundary {xd = 0} and 'purely elliptic' in the sense that it depends only on 
components 9m which are the arguments of elliptic profiles and is consequently elliptically polarized; these 
conditions on the error are described precisely with the hypotheses of Proposition 12 .401 Essentially, we form 
new solvability conditions which are the iterated profile equations modified by including the error terms. 

The second issue is that, since we do not make a small-divisor assumption, we can only guarantee 
solvability if we are working with finite trigonometric polvnomialj^. as opposed to general elements of H^^ 
with infinitely many nonzero Fourier coefficients. To resolve this, we approximate each vO^"+i by a finite 
trigonometric polynomial Vp'"^^, construct a corresponding corrector Vp, and use Proposition 12.391 and 
Lemma [2.421 to work with finite sums Up'"^^^ and Up'"^^ approximating U'^''"-+^ and Using these 

approximations, we are finally able to conclude (ll.63p . 



where I, O, V , andAf are as defined in the comments following Definition l 1.121 We also have the equivalent 
formulation Zi — {n £ Z : nluiio^ > 0}. 
(a) We also define 



2 Profile equations: formulation with periodic profiles 

Definition 2.1. (i) We define Z^' C Z^' by 

(2.1) Z^' := {a = {a,)fU ■ «^ S Z,}, 



where Zi is defined by 



(2.2) 




(2.3) 



{a £ Z : at most k components of a are nonzero}. 



Definition 2.2. For k = 1, 2, we define the following spaces: 



(2.4) 




For s > {d + 1 + 2)/2, it is clear that multiplication defines a continuous map 



(2.5) 



i?^^i(r|+' X X i/^'i(l++' X T*^) ^ H'-^rI^' X T^'). 



Remark 2.3. We define 



(2.6) 



C 



Ci X C2 X • • • X Cm, 



^'^For details on solvability, see Proposition 12. 91 
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where Cj is defined by 

( Rior i elUO, 

(2.7) C, := hlmz> 0} for i € V, 

[ {Im z < 0} for i e TV. 

For V e X T*f) where s > | + 3, one can show that spec V C Z^^'^ implies V extends holomor- 

phically in 6 to the interior of C*^. In particular, this uses the fact that then Ini(ai6'i) > for i ^ V U J\f. 
This allows us to make sense of 

(2.8) Uix,eo,^d) ■.^V{x,9o+uiiU...,9o+uiM^d) 
for 9o e T, > 0. 

2.1 The periodic ansatz and its corrector. 

Our periodic ansatz will have the form 

M Mm 

(2.9) v°(x, 9) = v{x) + 

m=l k=l 

a function in iJ'''^(M^^ x T^^), where s is to be specified, and which is holomorphic in 6* € C^^ C C*^, in 
particular. Using the notation 

(2.10) 9{9o,U) = {9o + WiCd, ...,9o+ uiM^d), 
a calculation shows 

(2.11) Kix) :=Z^°(x,0o,^<i)|,„^M,^,=ii = V°(a;,0(0o,Cd))|,„=^,5^^f^ 

results in the vanishing of the terms of order ^ when plugging into P{eu^, dx)ue in (ll.2l) (a). If we plug in 
a corrected approximate solution 

(2.12) 

<(x) {U''{x,9o,Cd) + eU'ix,9o,U)) = {V''ix,9i9oAd)) + eV\x,9{9o,U))) 

where we do not necessarily have of the form (12. 9p . but require V^{x,9) € H''''^{x,9), then the terms of 
order cancel out if and only if 

(2.13) £'{deo,d^,)U^ + L{dx)U" + M'{U")dg„U" = F{0)U°, 
using the notation 

(2.14) C'{dg,,d^J:=A{l3)de,+d^,, 

d-l 

(2.15) X'(Z^) :=^a„ij(0)Z^^j. 
A sufficient condition for (I2.13P is 

(2.16) Cide)V^ + L{dx)V" + X(V")9eV" = F(0)V°, 
where 

M 



(2.17) cidg) = J2 



m— 1 
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and 

M d-1 

(2.18) M{V)dg := 

m— 1 j—0 

Observe that the operator C{dg) is singular, so that one cannot simply invert it to solve for in p.l6p . 
One approach is to ensure the existence of a solution by imposing the following condition on V*': 

(2.19) L{d^)V^ + 7V/((V")5eV° - F(0)V" e Im C{dg), 

a condition which turns out to be equivalent to a differential equatiorF^ in V". 

The condition (|2.19l) on V" does not clearly determine y_ and the iJm.fe, the parts of V" which remain 
to be chosen. Strictly speaking, we will not satisfy (|2.16p or (|2.19p . Instead, we replace the differential 
equation (equivalent to (|2.19l) ') with another on^ll which agrees when we plug in = 9{0Q,(_ii) and can be 
decomposed into the system of profile equations (in v and the (Jm,k-) This allows us to solve for the profiles 
and a corrector satisfying an equation which agrees with (12.161) on 6* = 0(6*0, Cd)- Thus, this condition is 
also sufficient for (|2.13l) . 

Now we work towards defining the projection operators to appear in these differential equations, E and 
E'', which will allow us to construct the ansatz V" and its corrector V^. 

Definition 2.4. Setting <j) :— ((/)i, . . . , 4>m), we call a G Z*^'^ a characteristic mode provided det L{d{a-<j))) = 
and write a € C We decompose 

(2.20) C = u'^l^iCm: where Cm = {a G z'^'^''^ : a ■ (j) — Uaipm. for some Ua G Z}. 

Remark 2.5. Normally, the above is defined with Z^^'^ in place of Z^^'^, but in preventing our solution 
from blowing up, we have restricted its spectrum such that only a G Z^^'^ are considered. 

Now we are ready to define the projectors which will give us the profile equations and the differential 
equation for the condition p.l9|) . 

Definition 2.6. Let V G H^^'^ . The action o/E on V is defined by 

M 

(2.21) E = Eo + ^ E™, where EoV = and E„V = ^ P™T4(x)e"°^'", 
where Pm denotes the projection onto Ker L(d(f>jn), the action of¥^^ on V is defined by 

M 

(2.22) E^ = E^ + ^ E^^, where E^V = Vq and EJ„V = ^ PmVa{x)e'"-\ 

™=i QeCm\o 

and we use the notation 

(2.23) Eft=Eo+ H 

melUO meVUAf 

Remark 2.7. (i) It is shown in [2] that E is a continuous map in the H'' norm. Here, to complete the 
statement E : iJ'*'^ -> iJ'*'^, we must also have for each a G Cm that Ua G Zm- This is easily verified from 
the definitions of Zm, Cm, and the ri^. A similar calculation is done in Remark 12.121 (ii) It is not hard to 
show the continuity of E'' : TJ"-^ — > iJ"*-^. 

Remark 2.8. Denoting evaluation at 6 = 9{6o,^d) by $, observe $ o E = $ o E^ The projector E serves 
as the tool for solving for v{x) and the periodic profiles am,kix,9m) of our ansatz. The projector E^ is key 
to describing solvability (such as the condition (|2.19p ) and is thus used in the error analysis in solving away 
an 0(1) error output of $ o (J — E) written in the form $ o (/ — E''), so that we can prove Theorem 12.411 
That is the purpose which the following proposition serves. 

^^This differential equation is 1 12. 351 1. 
^■^This is the differential equation 1 12. 371 1. 



16 



Proposition 2.9. Given H G H^^(x,0) with finitely many nonzero Fourier coefficients, suppose 

(2.24) E^n = 0, 
Then there exists V e H^^{x,0) such that 

(2.25) C{de)V = U. 
Proof. We write out the series of T-L as 

(2.26) nix,e)^ H^{x)e'°''^ . 

Now we search for Va{x) such that, upon defining 

(2.27) V{x,e):^ J2 

one has 

M 

(2.28) jC{de)V{x,6)=J2 E ic^mL{dcj)^)Va,{x)e'"-' = ^ ff„(x)e^"-^ 
This holds if and only if for all a G Z*^-^ 
(2.29) 



where ut — (wj^, . . . , w^,/). We handle first the case where a e Cj \ for some j = 1, 2, . . . , M. Note that 

(2.30) 0^E)H{x,9)^ PjHc,{x)e'°'-\ 

aeCj\0 

which implies for each a G Cj 

(2.31) = PjHa{x). 

Recalling from Lemma Fl . 141 that Im A'^^L{d(j)j) = Ker Pj, we see there exists Wa{x) £ H^{x) such that 

(2.32) Haix) = iL{d(l3j)Wa,ix). 
Observe then 

(2.33) Ho^ix) = ^Z(n„/3,r^„w,) — M^a(2:), 



(2.34) = zL ( ^ a„/3,a-L^J ^ 

\m=l / ° 



Wc^ix). 



So we may satisfy (|2.29p by defining Kt(x) = ^M^ct(x). 

For the case a = 0, where EqH = iJa = 0, or for any other a such that Ha = 0, we may satisfy (|2.29p by 
taking Va = 0. 

For a ^ C with Ha 7^ 0, we have det L{J2m ^mf3, a ■ uj) ^ hy definition of C, so then p.29p can be solved 
directly for Va{x). Since there are only finitely many nonzero Ha, we find the same holds for the Va, and so 
the sum in ()2.27p indeed describes an element of H^^{x,6). □ 
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Remark 2.10. (i) In view of Proposition 12.91 in order to construct a corrector appropriate for our 
periodic ansatz V'' to satisfy p.l6|) . it is tempting to require something hke 

(2.35) E^(i(a^)V" +X(V")9eV°) = E^(i^(0)V°). 

However, note that to use Proposition l2.9[ if the we seek has infinitely many nonzero Fourier coefficients, 
(|2.35p alone is insufficient; we must use finite trigonometric polynomial approximations. Furthermore, even 
if V*' is replaced with a finite trigonometric polynomial approximation, there is no clear way to solve (|2.35p 
for V° of the form (|2.9I) . We explain: the components of the projector ¥^ are the E^; typically, one obtains 
some profile, say, <Jm,k{x,0„i) in (j2.9l) by examining the mth component of a projector. However, E^ maps 
into H^''^(x, 9), meaning its output generally varies with more than one component of 9 as opposed to varying 
with only the 6',„ component. Meanwhile, E has components E^ mapping into H'^{x,9,n), spaces suited to 
<^m,k{x,9m)- In fact, V° having the form (12. 9p is equivalent to having 

(2.36) EV" = V". 

(ii) With the considerations made in Remark |2.8[ it is natural to require instead of (|2.35p that 

(2.37) E{L{d^)V" +MiV^)deV°) = E(i^(0)V°). 

In fact, (j2.37p gives the solvability conditions discussed in Remark 11.151 which ensure the existence of a 
corrected approximate solution u'^{x) resulting in the vanishing of the terms of order 0(1) in (jl.2p fa'). To 
see this, recall that the 0(1) terms vanish if and only if (|2.13p holds. If we satisfy (|2.37p . we have 

(2.38) (/ - E)(i(9^)V" +MiV°)deV" - F{0)V°) = L{d^)V" + MiV°)dgV" - F{0)V°. 
Once V° is determined by imposing (I2.37|) . with Proposition 12.91 one can solve for a corrector in 

(2.39) Cide)V^ = -{I - E^)(Z(a^)V" + MiV°)dgV" - F{0)V°), 

(if we assume for the moment that we are working with finite trigonometric polynomials.) Although such 
does not necessarily satisfy (|2.16l) . after evaluating at 9 = 9{6Q,^d), which, recall, is denoted by $, it follows 
from (log)) and (pfg^]) that 

(2.40) $ mde)V^) = -$ (^L{d^)V° + MiV°)deV" - F(0)V") , 

where we have used the property $ o (/ - E'') = $ o (/ - E). In fact, for = <i>{V^) and U" = $(V"), ([^0)1 
is equivalent to (|2.13p . giving Ug(x) solving (|1.2l) (a) to order 0(1). 

It is from the components of (j2.37p that we will obtain the components cfm,k, ttl = 1, . . . k ~ 

1, . . . ,/i„ij of V°. Now we collect the set of equations which determines V°(x,6') € i/'''^(]R.'^^^ x T*^) for s 
sufficiently large (specified later): 

a) EV° = V° 

^2 ^) E (i(a^)V° + M{V°)deV°^ = E(F(0)V°) in > 

c) BiO)V%x',0,9o,...,eo) = Gix',eo) 

d) V° = in i < 0. 

In summary, (|2.41l) (a) gives that £(9e)V° = 0, and as a result, in the system for the original solution 
as a perturbation, this guarantees that (|1.2p (a') is solved to order O(i-). Equation (I2.4ip fb) represents 
the solvability conditions discussed in Remark 11.151 guaranteeing the existence of a corrected approximate 
solution Ug(x) as seen in (|1.3ip which solves (ll.2l) (a) to order 0(1). Finally, equations (|2.41l) (c) and (|2.41l) (d) 
straightforwardly correspond to (|1.2p (b) and (ll.2l) (c). 

Resonances. 

The next major task is to rephrase (|2.4ip (b') in terms of v and the am.,k to get the interior profile 
equations. (|2.4ip fb) will be broken down into components where E is replaced by Eq and Ep, p = 1, . . . , M, 
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yielding equations for v and the <Jp^i, I — 1, . . . , /ip. The nonhnearity in (|2.4ip fb) leads to resonances, which 
are generated in products such as aq,k{x, ■^)de,.cFr,k'{x, ^) when there exists a relation of the form 

(2.42) np4>p = nq(j)q + nr(j)r where p G {!,..., A/} \ {q, r} and (p, r) e Z x x Z,-. 

This means that (pq oscillations interact with (j)r relations to produce (j)p oscillations. While such <j)p oscil- 
lations arise on their own in products such as aq^k{x, ^)dg^(Tr.k' {x, not even some 9p dependence is 

present in the product aq^k{x, 0q)d0^ar,k'{x, Or) before we evaluate 9p — Oq = ^ , Or = Meanwhile Ep 
accounts for these interactions when it maps such {Oq, Or) oscillations to dp oscillations, as it maps 6) 
to H''{x,Op). When we express E_^ (|2.4ip fb') solely in terms of functions of {x,dp), it takes some effort to 
evaluate Ep{Ai{V^)dgV'^) in terms of the profiles, in particular to describe these interactions in terms of aq^k 
and CTr,fc'. This is the purpose of the interaction integrals, (I2.57p . appearing in the profile equations. To 
handle these difficulties we introduce some new definitions regarding resonances and interaction integrals. 

Definition 2.11. We say {(j)p, (j)q, (j)r) is an ordered triple of resonant phases and that {(f)q,(j)r) forms a 
p-resonance if 

(2.43) «p0p — nq(f>q + nr(l)r 

for some {np,nq,nr) Z x Zq x Zr, each entry nonzero. The triple {4>p, (l>q., 4>r) is called normalized if we 
have both (i) gcd{np,nq,nr) = 1 and (ii) ifpCzTUO, then rip > 0. 

First, we explicitly treat the case that the only normalized triples of resonant phases are permutations of 
{(j)p,4)q^(j)r) corresponding to one of the six rearrangements of (12.431) . It follows from the next remark that 
the sign of Up is determined, so this normalization uniquely determines rip, nq, rir- 

Remark 2.12. For a triple of resonant phases (0p, (f)q, (pr) with (j2.43p . one has {up, riq, rir) E Zp x Zq x Zr- 
To see this, taking the imaginary part of np(j)p — UqCpq + UrCpr, observe 

(2.44) nplm lo^ — Uqlm + rt,.Im w,.. 

Thus, if (0g,0r) forms a p-resonance, since ?iqlm and n^Im are nonnegative, so is rip Im oj^, so it 
follows that Tip G Zp. 

We now make a related observation with the following proposition which will greatly simplify our solution 
of the interior equations. 

Proposition 2.13. (i) A hyperbolic resonance (i.e. a p-resonance, where p E T U O) can only be formed by 
a pair of hyperbolic phases, and (ii) an elliptic resonance can only be formed by a pair including at least one 
elliptic phase. 

Proof, (i) Suppose for some p G TUO that np(j)p — nq(j)q + nr4ir with (rip, Uq, Ur) G ZpXZqX Zr- Considering 
the imaginary part gives 

(2.45) = riqlni + n,.Im w^, 

so recalling ni G Zi, rijlm uj_^ > 0, we see Im uj^ = Im cj^ = 0. Hence q,r G I U C The proof of (ii) is 
similar. □ 

2.2 The large system for individual profiles. 

Recall the decomposition of the projector 

M 



(2.46) E = Eo + ^ E„ 
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For m — 1, . . . , M, we enumerate by {£m,k, k — 1, . . . , z/„j} a basis of vectors for the left eigenspace of 



(2.47) tAW) = A^\0){rI + Y.MOh^) 
associated to the eigenvalue —i^„n and with the following property: 

(2.48) lm,k ■ rm',k- 



1, if TO = m' and k = k' 
0, otherwise 



For w e set 

(2.49) P™,fcV 
So now we may write 

(2.50) 

where 

(2.51) E„,fc(T4e 



v)'i'm,k (no complex conjugation here) . 



E = Eo+ ^ Y,^ra,k, 

m—1 k—1 



i^.e^ f(P™,fcK.)e"°^", a e C„ \ 
I 0, otherwise 



Recall the potential use of the projectors discussed in Remark 12.101 We will obtain a system of equations 
for v{x), am,k{x,Om) by applying Eq and the Em,k to (|2.4ip (b). 

We omit proof of the following lemma, which is almost identical to that of Lemma 2.11 of [2]. 

Lemma 2.14. Suppose EV° = V". Then 

(2.52) E™,fc(Z(9,)V°) - {X^ 

^ O'm ^k^T^m.k 

where X^^ is the characteristic vector field associated to (j)^: 



(2.53) 



We proceed by defining the interaction integrals which will appear in the interior equations. 
Definition 2.15. Suppose {(l)q,4>r) forms a normalized p-resonance, with 

(2.54) f^p'^^>p ~ nq(j)q + nr<j)r- 

For any f € H^{x,9q) define fn^ £ H^{x,9q) to be the image of f under the preparation map 



(2.55) 



f{x,0q) ^Y.M'^y' ^T^M^y'-'' 



kei 



Suppose s > + 1 and that cTq^k, <^r,k' G H^(R^ x T). We get exactly two normalized p-resonance 
formations from the arrangements of p.54[) ; 

(2.56) np(l>p = nq<pq + nr4>r, np(j>p — nr4>r + nq(j)q. 

To these equations we associate, respectively, the two families of prepared integrals; 
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(2.57) 



TT- / (cr<J,fe)n, X, —Op -Or de.CTr.k' {x, Or)dOr, fc G {1, . . . , fig}, fc' G {1, . . . , fir}, 



Jp^n,.,ng (x. Op) 



TT- I {crr,k)n^ \X,—Op -Oq ] de^CTq^k' {x , Oq)dOq, k ^ {1, . . . , ^r} , k' E {1, . . . , flq} . 

Zn ,/n \ rir TLr J 



Remark 2.16. Strictly speaking, Jp;,^^,„,, and Jp£,„^ are symbolically /„";_„p_„, and /„; ,„„^,„^ (as defined 
in [5], with respect to the triple {(f>q,4)p,(f>r) and nq(j)q = ripifip + n'^4>r) where n'^ — —rir. The difference is 
that we do not require nq>0, q<p<r- recall, instead, our concern is that (rip, Uq, Ur) & Zp x Zq x Zr- 



The following proposition shows that the prepared integrals 'pick out' the p-resonances. 
Proposition 2.17. Suppose s > + 1 and that Oq^k, cTr.k' G H^(R'^^ x T) have Fourier series 



(2.58) (Tq^k{x,Oq)= aj{x)e'^'^'' and ar^k'{x,Or) ^ J2 bjix)e'^'^-. 

jeZg\o ]ez,.\o 

Given {4'qT4'r) forms a normalized p-resonance, the prepared integral Jpn^ „^(x, 0p) belongs to H^~^{x,Op) 
and has Fourier series 

(2.59) Jp:nlnM,^) = E«J».(^)^^»^(^)* ■ (Jr^r)e''"'''. 

For the other integral in (I2.57p . one switches Uq and Ur above. Moreover, the junctions Jpn^ j^^{x,Op) can 
be described for Op G Cp through analytic extension of the expansions into the complex half plane Cp in the 
same manner that the profiles <Jm.k{x,Om) are extended into C,n in Section \1.2.1i 

Remark 2.18. We note that for p e V U Af, in (|2.59p . we can sum over just j G Z+ \ 0. To see this, 
observe that then (0^, 0^) forms an elliptic resonance, so Proposition 12 . 13l guarantees one of (t)q, 4>r is elliptic 
~ without loss of generality, say (j)q. Thus, Zq\0 consists only of negative integers or positive integers, and 
the only nonzero Fourier coefficients of (Jq^k are those aj{x) appearing in (|2.58p . with j £ Zq\ 0; we have 
aj{x) = for all other j. For example, if q €V, then Zq = Z+, and so, noting Uq > 0, we have the aj„^(a;) 
appearing in (j2.59l) are only nonzero for j > 0. Indeed, q Cz M implies the same of the aj„^(a:), so we may 
sum over just j e Z+ \ in cither case. We add that, for such j, the e*^""^" factors are bounded as Op ranges 
over Cp. 

Proof of Proposition \2.17\ The main part of the proof is showing that the expansion (|2.59|) converges to the 
desired result in H^~^ (x,Op), where we have restricted our attention to real Op. The proof of this is very 
similar to the proof of Proposition 2.13 from [5]. Analytic extension into the complex half plane Cp is done 
in the same way that the expansions for the profiles are extended into their corresponding complex half 
planes, as discussed in Section ri.2.11 □ 

Interior equations. 

By applying Eq and the Em,k to (|2.4ip fb). we obtain the system for v{x) and the am,k{x,Om) with the 
following proposition. 

Proposition 2.19. Suppose {(j)q,(j)r) forms a normalized p-resonance, with 

(2.60) rip(f)p = UqCfiq + nr4>r. 

and that all other normalized triples of resonant phases are permutations of {(f>p, (j)q, (f>r) each corresponding 
to one of the six rearrangements of (|2.60p . 
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(i) Given V" as in (j2.36l) and represented by (j2.9p . the equation (|2.4ip fb) is equivalent to the following 
system: 



(2.61) mM^Y^Y. E ^ / 

j=0 m=l k,k' = l ^"^0 



(2.62) 



(a) X0^CTp,/(a;,6lp) + XI E °^ii,M)^Ov^P^k'{x,ep)+ 

d-l P-p d-1 Mp Mp 



Li— -L A-p ^-p U— i A-p A-p -j^ />Z7r 

EE E ^p,tj'^p.'=(2^'^p)^«P^p,fe'(2^'^p) "EE E ^p>02^ / <^pA^^^p)9ep<^p,k'{x,ep)dep 

j=0 k=lk' = l j=Ok=lk' = l 



Mp 



(c) Jp^i.j (x, 6'p) = X (CTp^fe [x^Op) 



fc=i 



w/iere 



(2.63) 



M£ p.r Pr_ Pq 

J'M jk,k' a -\ \^ 



Jp,Lj{x,9p) = < 



fc=ifc'=i fe=ife'=i 

Alp Mr /^r Mp 

E E '^q:-n,,nAx^Sp) + E E ^pL Jg:t,-nJ^'^p)> P^g' 
fe=lfc' = l fe=lA:' = l 

Pq^ Pp_ Pp_ Pq_ 

E E '^L,nq,-nj,{x,Op) +X X '^r^,-nj,,ng{x,9p), P — Z, 



k=l fe' = l 

0, otherwise; 



k=l k' = l 



(2.64) 



where 



(2.65) J^p,/,j(x,0p) 



d-l Mp 

(a) X0^crp,,(a;,6'p) + ^ ^ ap';,j(ll)<9epCrp,fe'(a;, 6*^) + 

j=0 fc' = l 

d-l Mp Mp 

(^)EEE 

j=0 k=l k' = l 

d-1 Pp 

(c) Y Jp^l,j (a^' ^p) = E Sp,/^P.fe ^p) 

i=0 k=l 

forpeVyjN, le {i,...,fip}, 



{Pq Pr Pr Pq_ 

E E s£ 4:^,nM ^p) + £ £ ^,C,n, (-^ ^p)' p = 
0, otherwise. 



k=l k' = l 



The constant vectors -R*^',^j appearing in (j2.6ip are given by 



(2.66) 



R'jjn = /3j(9„Aj(0) • r„^^k)rm,k'- 



The constant scalars bp'ij, ^y'lj^ '^p'lj' ^pi' '^^^ coefficients of the scalar linear function of v, a'^ij{v), 
in (|2.62p - (l2.65p are given by similar formulas, but now involving dot products with the vector £p^i. 
(ii) Equations \2.Q1\ . ()2.62p . form a hyperbolic subsystem, that is, a system in v and CTp^i for p G X U O 
independent of Uq^k for all q V U Af. 
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Remark 2.20. While the hyperbohc interior equations are independent of the elhptic profiles, the elliptic 
interior equations (j2.64[) are not in general independent of the hyperbolic profiles, since these appear in the 
elliptic interaction integrals if some phase paired with a hyperbolic phase forms an elliptic resonance. 

Proof of Proposition \2.19\ Regarding (i) , (12.611) follows from application of Eq to (|2.4ip (b) , and (|2.62p from 
application of Ep_/ to (|2.4ip fb). One similarly arrives at (|2.64p : in doing so, one finds that the second triple 
sum appearing in (12.621) fb) is zero for p S T' U M , since then spec CTp.fc, spec Cp^fc' C \ which is either 
Z"*" \ or Z~ \ 0, resulting in (|2.64l) (b'). The differences between (|2.63p and (|2.65l) merely account for the fact 
that for p G I U O, we have rtpGZp\0 = Z\ implies —rip e Zp \ = Z \ 0, while for p e T-" U A/", we have 
rip £ Zp\Q implies —rip ^ Zp\Q. The reverse implication of (i) follows upon recalling the decomposition 
(12301) . 

For proof of (ii), first we show the elliptic profiles drop out of (|2.6ip . This follows from the argument made 
in proving (i) that terms such as those in the triple sum of (|2.6ip with m G V U Af are zero, since then 
spec (Tm.k, spec crm,k' G Zp\0. Finally, to see that no <7q^k with q £ V yj M appears in (|2.62p . we have only 
to check the interaction integral terms Jp.i^j of (|2.62p ('c). However, the terms of (|2.63l) defining Jp,;j involve 
only profiles corresponding to p-resonant phases; Jp,/j- is only nonzero for some p G I U O if p = p, or r 
and, by Proposition 12 . 1 31 p, q, r G XU C Thus, in this case, no elliptic profiles appear in Jp,ij- □ 

From the above, we find the hyperbolic subsystem has the same form as the system in [5] and we will 
solve it similarly. On the other hand, we will not solve the elliptic interior equations exactly. 
Boundary equations. 

Now that we have the interior equations for the profiles, we formulate the boundary conditions to be 
imposed on the profiles. From (|2.41l) (c). we get 



(2.67) 



(a) B{{))v = Gix') 
(6)B(0)V°*(x',0,^?o,...,^?o)- 

^ ^o,n,k(x',0,ea)r„„k+ ^(7™,fc(a;',0, 

meXuVUAf k=l meO k=l / 



where V°* denotes the mean zero part of the periodic function V°, and we have similar for G* . While ()2.67p (a) 
gives boundary data for v, it remains to establish satisfactory boundary conditions on the individual profiles 
Cm.fc from (|2.67p (bV This is done with Proposition 12.221 

Lemma 2.21. Each of the sets of vectors 

(2.68) {B(0)r„.,fc : m e I U V , k = I, . . . , ^i„,}, 



(2.69) {B(0)r™,fc : m G X U AA, fc = 1, . . . , /x^}, 
is a basis of . 

Proof. First we prove (I2.68P is a basis of C^. Recall our assumption of uniform stability. Assumption (11.51) . 
which tells us that B{0) maps a basis for the stable subspace E'*(t, r/) to a basis for C^. Thus, according to 
Lemma ri.l3[ which states that 

(2.70) E^(t,^) = ©™eiupKer L(d0™), 

we have that the set (|2.68p is a basis for C^. To see the same holds for (|2.69p . recall from Remark If. lOf ii) 
the fact that for each nonreal lo^, say, without loss of generality, m d V, there is another eigenvalue w^/ 
satisfying = thus with m' G Af, and that similar holds for the associated eigenvectors. It follows 
that 

(2.71) {S(0)r„,fc:mGlUAA, /c = l,...,/x„} = {S(0)r™,fe : m G X U 7^, fc = 1, . . . , 

(2.72) = {S(0)r„,fc : m G I U 7^, fc = 1, . . . , fi^}. 

Clearly, p.72p spans the same subspace as p.68p . i.e. C^, and so (|2.69p is also a basis of C^'. □ 
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Proposition 2.22. The data {(Jm.k(x' , 0, 6*0); m E XUV U Af, k G {1, . . . , /i,n}) are determined by the data 
(cm,fc(2;', 0, f?o); E 0,k E {1, . . . , /i„i}) in that there exist constant matrices such that the zero-mean 
boundary condition, 

(2.73) B{0)V"*{x', 0,00,..., eo) - B{0) ( E '^-^-k-A^' ' &o)r„^.k \ = G*{x', 0o), 

\m=l k=l / 

is equivalent to the condition 
(2.74) 

i'^i.ki^'^ 0' ^0); mElUVuM,kE{l,..., fim}) = M± iG*'^,a^^Jx', 0, Oo); m E 0,k E {1, . . . , /x„}), 
(using + (— ) to denote a part with positive (negative) spectrum.) 
Proof. We rewrite (|2.73p as 

B{0)V°*{x\0,9o,...,9o) = 

(2 75) / \ 

Oo)rm,k — G*{x' , 9o). 

\7neXUVUAf k=l meO k=l J 

Recall "P-profiles have positive spectra and A/'-profiles have negative spectra. Thus, using the subscript n to 
denote the nth term in a Fourier series, we get for n > 0, 



(2.76) 5(0) ( E'^'".'^-.»(^''0'^o)^'".'^- I = Glix'M - ^(0) ( ^ E 

^m^k,n(^X , 0, ^o)^m,/i: 

and for n < 0, 



\meIU'P k=l J \meO k=l 



(2.77) B(0)[ J2 E^".'=."(^'''^'^o)'^™.fc) =^"(^''^o)-^(0) [ E E^™'fc'»(^''0'^^^^ 

\meIU7^^fe=l / VmGO fc=l / 

By LemmalMIl both sets {B{0)r,n,k : fc e {1, . . . , /Xm}, m E lUV}, {S(0)r„,fc : fc e {1, . . . , ^i™}, m E luW} 
are bases for C^, and the desired result follows from this. □ 

So we now have a sub-system f (l2.6ip . (|2.62p ) in just the hyperbolic profiles am,k, m E lUO, and the mean 
V with boundary data for v and determination of X-boundary data from O-boundary data as in ()2.67p fb). 
where (fl!67l) (b) is equivalent to (f2J4|) . 

Remark 2.23. We will seek v, a-pj such that, in addition to the equations above, we satisfy the initial 
conditions 

(2.78) v = and apj = in t < for ah p, /. 

The large system consists of the interior equations (|2.6ip - (|2.65p with boundary equations (|2.67p and the 
initial conditions (|2.78p . 

2.3 Solution of the hyperbolic sub-system of the large system. 

First we will obtain the 'hyperbolic part' of our solution, V^, which will satisfy, in place of ()2.4ip . 

a) E„V,*J = V° 

c) B(O)V^(a;',O,0o,...,^o) = G(a;',0o) 

d) V^^Omt <0. 
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Condition (|2.79p (a) serves the same purpose as (|2.4ip fa'). but adds the restriction that V° only consists of 
the mean and the hyperbohc profiles. Thus, recalling the form for V°, we see V° takes the form 



M 



(2.80) vl{x,e^,...,eM) = v{x)+ E 

meXuO fc=l 

The condition (j2.79p fb') is exactly the hyperbolic sub-system f (|2.6ip . (|2.62p ) found in Proposition l2.19l This 
follows directly from Proposition 12.191 and the definition of E/j . Since V^jl has the form (|2.80[) , the boundary 
condition (|2.79p (c) becomes (|2.74p . ignoring the appearance of the elliptic profiles {<Tm,k with ni ^ V U Af) 
in the left hand side. 

To solve this, we employ the same approach as that which is used in p] to solve the full system of profile 
equations. This is done by solving the system with an iteration scheme such as the following 

a) E,Vr = 

b) Eh (i(a,)V,';^" +X(V°'""^)a,V,';^") = E,(^^(0)V°'"-^) in > 



(2.81) 



c) i?(O)V°'"(x',O,0o,...,eo) = G(a;',0o) 

d) V°'" = in i < 0, 

where we note that (|2.8ip (a) means V'^'"' is of the form 

(2.82) V,^'"(x,0l,...,^^M)=«"(a;)+ J2 E ^")^™.^' 

melUO k=l 

where iterates w"(a;), cr^ ^.(a;, 6m) for m £ X U O satisfy an iterated version of the hyperbolic sub-system of 
the profile equations, f (j2.6ip . (|2.62l) ) which is encapsulated by p.Sip fb). 

The following proposition gives the solution to the iterated system as well as the solution to the hyperbolic 
sub-system itself. We omit the proofs of these as they are almost identical to the proof of Proposition 2.19 
together with that of Proposition 2.21 from [2]. 

Proposition 2.24. Let T > 0, m > ^ + 1 and suppose that G(x',9o) G HJP". 

(i) Setting V^' = 0, there exist unique iterates V^'" G H^, n>\, solving the system (j2.8ip . 

(ii) For some < Tq < T the system (|2.79p has a unique solution G ^Tq- Furthermore, 

(2.83) lim V°'" = V? mM"'-\ 

n— >-oo 

and the traces V^\xa=Q) ll\xd=0! '^^'^ o'p,i\xa=o, p T U O, all lie in HJp^. 

2.4 Approximate solution of the equations for the eUiptic profiles. 

Recall from Proposition 12.191 the elliptic interior equations: 

rf- 1 Mp 

(a) X^^apj{x,ep) + X! E °^P,ijfe)"5epC^P,fe'(a;,6'p)+ 

j=0 k' = l 

d—l A^p Mp 

, {h) EE E ^p,tj'^p,fc(2^'^p)'9epf^p,fe'(2^,6'p)+ 

(2-«4) j=ofc=ife'=i 

d-l lip 

(c) ^ Jp,/ J (x, 6lp) = ^ ep_,crp,fe {x,ep) 

j=0 k=l 

iovpeVuM, I e {1,. . .,/ip}. 

These complex transport equations may not generally have exact solutions. Instead, our elliptic profiles dp,; 
will approximately solve these equations in the sense that they will hold at the boundary Xd = 0. It is 
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sufficient to simply evaluate the above expression at Xd — 0, isolate dxj^crpj{x' ,0,9p), and require that apj 
has the appropriate Xd-derivative at the boundary. Of course, we will also have to adhere to the boundary 
conditions, so that the trace ap^i{x' ,0^9o) at the boundary satisfies (|2.67l) (b). and we must satisfy initial 
conditions (j2.78p . With Proposition l2.30[ we construct such elliptic profiles ap^i{x,9p) by solving some wave 
equations in which Xd plays the role of the time variable and with appropriate boundary data corresponding 
to ap^i{x' ,0,9p) and d^^apjix' ,0,9p). 

In addition to the exact elliptic interior equations, we consider semilinear equations in iterates (7^; for 
p & V U A/", I S {1, . . . almost identical to the iteration scheme which is used to obtain the hyperbolic 

profiles. While we do not need elliptic iterates ct^ ^ to get the elliptic profiles (Jp^i, which are instead obtained 
by the process described above, they allow us to handle hyperbolic parts and elliptic parts uniformly in the 
error analysis. The semilinear equations in the elliptic iterates a^i are 

(2.85) X^^<j^^i{x,9p)+ 



(2.86) 



d-l lip 

EE 



- 1 f-^p l-ip 



j=0 k=l k' = l 



(2.87) 



EE2. 

j=a fe=i fe'=i 



Alp 



1 Mp fip 



"p,l,j P,nq,n. 



.(^.^p) + EEE 

j=0 k=l fc' = l 



jfc.fc' jk,k'. 



(X,9p) 



IJ.p 

(2.88) =E4'<fe'(^'^p)' forpeT'UAA, /£{!,..., /xp}, 

fe=l 

where we define 

(2.89) Jp:nSM dp) = ^ f\^';:k'u U ~ r^^O (^r)d9r, 

and similar for Jpn^'n^- Imposed on the iterates a^i are boundary conditions identical to those for the 
profiles (Tp^i, i.e. p.74p with a"^ j. in place of each CTm.fc, and initial conditions 

(2.90) cr;^^, (a;) = for t < 0. 

Similar to the elliptic profiles, the elliptic iterates (7^ ; will only approximately solve (|2.85p - p.88p in the sense 
that these equations will hold at = 0. In fact, the elliptic iterates a^^ are obtained with the same kind 
of construction used for the elliptic dp,;. In addition to the elliptic profiles, their corresponding iterates are 
constructed in Proposition 12.301 

For now, let us fix p e {1, . . . , M}, I G {1, . . . , /Xp} and an integer n and take the task of finding approx- 
imate solutions CTp^i and a^i of (j2.84p and the system (|2.85l) - (|2.88p . respectively. First we rearrange (|2.84l) . 
isolating the vector field applied to ap^i in the left hand side, and denoting the resulting right hand side by 
fp^i, getting something of the form 

(2.91) X^^ap^iix,9p) = fpd{x,9p). 

Similarly, we isolate X^pO-p ; in the left hand side of (|2.85p - (|2.88p . and define fpi{x,9p) to be what remains 
on the right hand side, so (I2.84p becomes 

(2.92) X^^a;,{x,9p)^f^,{x,9p). 

Now we isolate the quantities this prescribes for the traces of the C7p_; and the CTp ; at the boundary Xd = 0. 
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Definition 2.25. (i) The coefficient of dx^Cp,! in the left hand side of (|2.9ip is one, so we may isolate it in 
this expression. We will require this equation to hold at xa = 0, and define bp^i so that doing so is represented 
by the condition 

(2-93) dx^cTp,i\xi=o — bpj, 

noting that the right hand side depends only on the boundary data <Jm,k\xd=o-i m = 1, • ■ • , Af, fc = 1, . . . , 
(a) We define ap^i such that the boundary condition on Op^i in (|2.74p is equivalent to 

(2-94) (yp,i\xd=Q = a,p,u 

noting that the right hand side has already been determined by our solution of the hyperbolic profiles. This 
thus determines the right hand side of (|2.93p . 

(in) We isolate dx^f^^^i on the left hand side of (|2.85p - (|2.88p . and define 6^ ^ such that requiring this to 
hold at = is equivalent to 

(2.95) dxdali\xd=o^hli. 
(iv) Consider 
(2.96) 

0,d„);melUVuN,ke{l,..., m™}) = M±(G*^± , ^^(x', 0,d„);m e 0,k £ {I, . . . , ^„J), 

for as defined in the proof of Proposition 12.221 We define ^ such that the condition on ; in (12.96^ 
is equivalent to 

(2.97) <zU.=o = a;;. 

Lemma 2.26. Suppose G € . Let V° € ^t^^ be the solution of the hyperbolic sub-system constructed 
in Proposition \2.24\ For a^^, b^i, Gp^i, and bpj as defined in Deflnition \2.25\ we have a^i G H^^(x' ,9p), 
bp A e H^^{x',9p), and ' ' 

(2.98) lim a", = Opi m HS+\x',0p), 

(2.99) Mm b^i^bpj inHI.(x',9p). 

Proof. Checking a^i G H^^^{x\9p) is straightforward, since G{x',9p), ^(x', 0, ^p), m G O, are in 

H^+\x',ep). 

Now we show b^i G H^^{x' , Bp). Observe b'^i consists of the terms evaluated at = (g^, (g^, (g^, 
and 

(2.100) [dxd~X^^)cjli{x',Q,ep). 

Since the i;"^, cr" j,(x', 0, 0p), p G XU O, are in H^^^{x',ep), and H^^{x',0p) is a Banach algebra, it follows 
that (|2.100| 'and both ([2?86| and ([ZM]) at Xd = are in H^^{x' ,9p). It remains to show this for (|2?87| at 



Xd = 0. This follows from Proposition 12 . 1 71 with x' in place of x. Taking into account that, for p G ZU O, 



(2.101) lim <,.U,.o = ^p,fcU.=o in H^+\x',ep), 



(2.102) lim dxdCrp^k\xd=o = dxdCrp,k\xd=o in H^„{x',0p), 

the proof that ((2J8)) and ((2^ hold is similar. □ 

Remark 2.27. To simplify obtaining our approximate solutions of (|2.85p - (|2.88p . it will be useful to extend 
functions a^i^Up^i G H^^^ = H''+^{{-oo,Tq) x E'^-i x T), b'^i^bp^i G H^^ = H^{{~oo,Tq) x M'^^i x T) on 
the half-space to elements of H'+'^iR'^ x T) and i7"(M'^ x T), respectively. 
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Lemma 2.28. For s > there is a continuous extension map 

(2.103) E : H'{{~oo,To) x W'-'^ x T) ^ iJ-'(M'' x T). 

Proof. It is shown in 4.4 of [9 that there is a continuous extension map from i7'*(K'|_) to i?''(R'^). □ 

From now on, in place of ^ , Op^; , 6^ ^ , bpj we refer to their extensions to the respective spaces mentioned 
above unless we exphcitly state otherwise. 

The fohowing lemma is a standard kind of result in the theory of hyperbolic initial/boundary value problems. 

Lemma 2.29. (i) For a e H'+^iW^ x T) and b e H''{R'^ x T), there exists unique <r e H'^^ = H''+^{R'^ x 
[0, D] X T) solving the wave equation 

(2.104) a2^^-A,,,e„^ = 0, 
and initial-Xd conditions 

(2.105) cU,=o = a, 

(2.106) 5,.^U.=o - b. 

(ii) Furthermore, <; satisfies 

(2.107) kl^^+i <C(|a|ff.+i + |6|H0, 
/or some constant C independent of the choice of initial data {a, b}. 

Justification of Lemma [2.291 can be found in Chapter 6 (see 6.18 and Remark 6.21) of [3]. 

Proposition 2.30. For a^i, b^i, Op^i, and bp^i as defined in Definition \2.25[ where p Cz V U M there exist 
api{x,9p), ap^i{x,9p) € H^^{x,9p) with compact Xd-support in [0,1?] satisfying 

(2.108) a^i = <jp^i = 0, for t < 0, 
boundary conditions 

(2.109) crp,i\xd=a = o,p,h crp,i\xd=o ^ a,p,h 

dxd'^p^llxa^O = bpi, dxa'^p'l\xd=0 = bp,l, 

and 

(2.110) lim a; I = apj m H'^+Ux^Bp). 

Proof. We apply Lemma 12.291 to initial data {ap^i^bpi + dtOpi] to obtain the corresponding solution of 
(|2.104p . denoted (,p^i € , and subsequently to each {a^ijb^i + dta^i}, obtaining solutions <jp; S H^^. 
Now we fix a smooth cutoff function xi^d) supported in [0, D), and define 

(2.111) apj{t,x'\9p) = x{xd)^p,i{t- Xd,x",9p), 

(2.112) a;^i{t,x'\dp) = x{^dK.i{t-Xd,x\dp). 

It is easy to check that then a^^, apj also belong to H^^. Since <^p,i\xd=o, <ipi\xa=o are supported in t > 0, 
by finite speed of propagation for solutions to the wave equation p.l04p . for any fixed r > 0, the {t,y,9p)- 
supports of <ip,i\xa=r, <ipi\xa=r are contained in the union of balls of radius r about points in {{t, y, 9p) : t > 0}. 
Thus, the support is contained in {(t,y,9p) : t > — r}. It follows that crp,i{t,y,Xd,9p), i{t,y,Xd,9p) are 
zero for all t < 0, since the supports of dp^;, CTp ; consist only of points satisfying t — Xd > —Xd, i.e. t > 0. 
It is easy to check from (|2.11ip that the ap^i, (jpi have the desired traces at Xd — 0, satisfying p.l09p . To 
establish (|2.110p . first we note we may regard the Cp,/, cTpj as elements of H^^{x,9p) by defining them to 
be equal to zero for Xd > D and restricting to t < Tq. One easily verifies the bound 

(2.113) |o-p_, - <Tp^i\h!^+^ ^ l^p.i ~ ^P:'lif^+i < C'lkp,/ - 
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Applying Lemma [2.291 to the solutions ; — (^pj, from the above and the estimate (|2.107p we get 
(2.114) - < C2 {\al, - ap,,U=+^ + W^i - bp,i\H^) . 

Now we conclude (|2.110p from Lemma 12.261 and the continuity of the extension map from Lemma 12.281 □ 

Now that we have obtained the elliptic profiles, we have determined the ansatz V*^. 

Definition 2.31. Let G G and let V° G ^t^^ be the solution of the hyperbolic sub-system constructed 

in Proposition \2. 24\ and let the V^'" S H^^^ be the corresponding iterates. Additionally, let Cpj, i S ^Th^ ' 
where p €V UM, be the elliptic profiles and iterates obtained in Provosition \2. 3^ 
(i) We define the elliptic part of our ansatz 



M 

meVUAf fe=l 



M Mm 
meVUAf k=l 



,k- 



(2.115) V°(a;,0i,. 
and the corresponding nth iterate 

(2.116) v"/\x,ei, 

Our ansatz is defined to be 
(2.117) 

and we also define V"'" V,"'" + V°'". 

(a) We isolate the error in the iterated interior profile equations, defining R"^{x,0) by 

(2.118) i?" := E{L{d^)V"'" + M{V^'"'^)deV°'" - F(0)V"^"-i). 
Proposition 2.32. Suppose the hypotheses of Definition \2.3l\ are satisfied. Then V° G H 

(2.119) lim V°'" = V° m H^„. 

n— >-oo 

Proof. The claims follow directly from Proposition 12.241 and Proposition 12.301 
Remark 2.33. Recalling ((2?8T|) . note that 

(2.120) E,,(L(5^)V°'" + X(V°'"-i)9eV°'" - F(0)V°^"-i) = 0, 
and so i?" of (|2.118p is purely elliptic in the sense that Egi?" = R" . Moreover, 



, and 

J-o ' 



□ 



(2.121) 



pevuAf 1=1 



which is zero at the boundary Xd — 0, since we have satisfied (|2.9ip for Xd — 0. Summarizing the results of 
Proposition 12.241 and Proposition 12.301 we conclude 



(2.122) 



a) EV°'" = V°'" 

b) E (Z(5:,)V°'" + X(V°^"-i)aeV°'") = E(F(0)V°'"-i) + i?" in > 



c) i?(O)V°'"(.T',O,0o,.-.,eo) = G(a;',0o) 

d) V"'" = in t < 0. 

The sense in which an error such as i?" is small is clarified by Proposition 12.401 
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Remark 2.34. Recall the discussion in Remark ll.lOf ii) which gives the bijection between the eigenvalues 
a;,„ with m £ V and w^/ = m' G Af and between the corresponding eigenvectors. A careful look at 
the profile equations shows that the equations for the elliptic profiles come in conjugate pairs. That is, 
the equation for am,k, some m Q V, is the conjugate of the equation for am',k' with the corresponding 
m' € Af and k' . As a result the solutions we have obtained, the elliptic profiles, also come in conjugate pairs, 
satisfying <Jm',k' = <^m,k- Meanwhile, the hyperbolic part of the solution Vj] is real. It follows from these 
observations and the definition of our ansatz that when we plug in 6 — 0{Oo^^d)^ getting 

(2.123) U°{x,eo,Cd) ■^V%x,0o+^,U---,Oo+^mU), 
we have a function U^{x, 6o, ^d) which is in fact real on w"^^ x T x M+. 

2.5 The expansions for the approximate solution 

For this study, we define the spaces and £^ as in [5]. 

(2.124) E^ = C{xd,H^{x\9o))nL\xd,H^+\x',eo)), 

where by C{xd, H^{x' ,6o)) we actually refer to functions in C{xd, H^(x' ,9o)) with x^-support in [0,1?] for 
some large enough D, and for C{xd,H^{x' ,9o)) we use the L°° (xd, Hj.{x' ,9o)) norm where the supremum 
is taken over Xd > 0. These spaces are algebras and are contained in for s > ^ii. Theorem 7.1 of [TT], 
as discussed in Section 11.1.31 regarding existence of solutions of the singular system, tells us that for s large 
enough, one has existence of solutions to (|1.36p in the space E^ on a time interval [0,r] independent of 
e £ (0,eo]. For these reasons, the space E^ and related estimates are key to our analysis, in particular for 
our main theorem (Theorem 12. 41[) which relies on Proposition 12. 43[ also proved in [TT]. Additionally we use 
the spaces 

(2.125) £^ = {Uix,9o,U) ■■ sup\Ui-,-,U)\E^ < oo}, 

id>o 

which also play a role in Theorem 12.411 The proof of Theorem 12.411 uses estimates regarding functions such 
as V in H^'^^ and the corresponding U in £^ which results from the substitution 9 = 9{9o,S,d)- While this 
moves us out of the space of periodic profiles H^^{x,9), we get functions in £^{x,9Q,£^d) which we can 
approximate with finite trigonometric polynomials (truncated expansionj^) in 9 = 9{9o,^d)- 

Definition 2.35. For k ~ 1, 2 we define 

(2.126) £f := {U{x,9o,U) = V{x,9)\g=e(ea,U ■ ^ ^ H't^^*), 

with the norm | • |f = . (Note: The subscript T has the same indication as it did for the H!^^''^ used in [5], 
where the subscript T is at first ignored.) 

It is verified in Proposition 12 . 39l that elements of £^'^ are bounded in the £^ norm. 
Convergence of expansions in £^. 

The following notation gives us a way to sort the spectra a G z^'''^ of elements V in H^^''^, which will 
aid in showing U{x, 9q, ^d) = V{x, 9)\g^g^g^^^^-f has an expansion converging in £^, with Proposition 12 . 391 

Definition 2.36. Let enumerate the set {a • a; : a G Z^^'^^}. For each j, k, we define 

M 

(2.127) := {a G Z^''-^ : ^ a, = j, a • ^ = Zfe} 

i=l 

and pick out one element G Cj^k ■ 

expansions are given by Proposition 12 . 351 
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Remark 2.37. A nice consequence of the fact that we work with ^^^-2 instead of general Z^''''' is that we 
can easily show the Cj^k are finitely. To see this, first fix j e Z, /c S {1, 2, . . .}, and p, q € {1, . . . , M} and 
take arbitrary a € Cj^k such that all but the pth and qth components are zero (either of the pth and qth 
components may be zero, as well.) We claim this is the only such element of Cj_k- This is because 7^ co^ 
implies 

(2.128) P ^ )( ^ 

has a unique solution (ap,ag). Thus, the number of such a in Cj^k is bounded by the number of pairs of 
components, so |Cj,fc| < M{M ~ l)/2. 

A fmiction V{x,e) e i?^+^'^(x, 6*) of {x,e) e M'^+i x C^^ has a series 

(2.129) Vix,e)^ J2 ^«(^)e'"'^ 

and, for fixed Xd, squared H^{x',9) norm 

(2.130) \Vix,0)\hi^',e)^ E E + 
From Sobolev embedding and the fact that 

(2.131) H^+\x,0) C LHxd,H^'^\x\0))nH\xd,H^{x\0)), 

we find V{x,0) e L'^{xd,H^'^^{x' ,0)) n C{xd,H^{x' ,0)), implying the partial sums of the series (I2.129P are 
bounded and converge in H^{x' ,0) uniformly with respect to Xd > 0. We will prove Proposition I2.3i?l bv 
using these facts with the following lemma, which shows for finite truncations of expansions (|2.129p that the 
norm | • \H^(x',e) dominates He=e(eo,5d)|^/j,(a.,.(,„) independent of {xd,£,d)- 

Lemma 2.38. Suppose V G H^^'"^ {x,0) with series given by (j2.129p . For 0{0o,£,d) o.s in (I2.10p and integers 
Ml < M2, Q < Ni < N2, we have the following inequality: 

Ah N2 M2 N2 

(2-132) IEEE K.(.)e-«(^-«''^|?,,(.,,„)<| E E E ^"(-K-%,(.',.)- 

j=Mi k=Ni aSCj.fc j=Mi k=Ni aGCj.fe 

Proof. We estimate 

M2 N2 M2 N2 

IEEE y<^(-y''-'''-'^%,i.',o) - lE E(E yo^i-y-'^y'''\hi^',ov 

j=Mi k=Ni a<£Cj,k j=Ah k=Ni aeCj,k 

M2 N2 

(2.133) = EEkE E ^f'^"(-)^'"'^'^)iW')(i + iJi)'^^"'"^ 

j = Mi \p\<s k=Ni aeCj,k 
Ah N2 

(2.134) < E E E E i^f'^"(^)^'"-"'''i'^(.')(i + i^i)'^^"""- 

j = Mi k=Ni aeCj,k |/3|<s 

For (|2.133p . we used the formula which that in ()2.130p generalizes. For a € Cj^k, we have |j| = | J2i — \^\ 
and Im(Q; ■ lu) > 0, so for all ^d ^ 0, the sum in (|2.134p is bounded by 

Ah N2 Ah N2 

(2-135) IEEE yo^i-y^-'\hi.;e)=Y. E E E l^f'^"(-)li^M(l+l"l)''^~'"^ 

j=Mi_ k=Ni af^Cj,k j=Ah k=Ni af^Cj,k |/3|<s 

where the equality (|2.135p follows from the formula in p.l30p . □ 



^*We work with Z^^'^ due to the fact that we only have to consider quadratic interactions. Further discussion on interactions 
and resonances can be found in Sections 12.11 and 12.21 
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With the estimate of Lemma I2.38[ we are ready to prove that elements of H^^'"^ yield elements with 
expansions converging in £^ upon the substitution = 0{Oo,^d)- 

Proposition 2.39. Let V G H^^''^{x,9) with expansion ()2.129p and set VI = V\e=g(^gg^^^j for 0{9o,^d) as in 
(PJU)) . Then we have U € 

(2.136) iZ^lfj <C|V|^^+i, 

and we have convergence in £^ to lA of finite partial sums, independent of arrangement, of the series 

Proof of Proposition \2. 39\ We explicitly show the convergence of the finite partial sums 
(2.138) ^ K(x)e*"-^(''°'«<^). 

— n<j<n a^Cj k 
l<k<n 

For integers ni, n2 with ni < "-2, an application of Lemma 12.381 gives 



(2.137) U{x,0o,U)^ J2 ^"(^'e" 



(2.139) I Yl E V»ix)e''^-'^'''^^^%^i.',Oo)<\ E E ^"(^'^ 



"l<b'l<"2 aeCj^k ni<b'|<"2 aeCj^k 

ni<k<n2 ni<k<n2 



It follows from (|2.139p that since the sequence of the partial sums 



(2.140) J2 E 



l<k<n 



is Cauchy in H^{x',0) uniformly with respect to {xd,£,d), so is the sequence (|2.138l) in H^{x' ,0q), and thus 
we have convergence. We similarly get convergence of the (j2.138l) in L^(xd, H^^{x' , 0o)) after integrating the 
inequality (|2.139l) (with s + 1 in place of s) with respect to Xd and noting the right hand side is independent of 
^d- Hence, the (|2.138p converge in It is not hard to show the limit is in fact V\g=g(^gg^^^'), and the estimate 
(|2.136p easily follows. The same proof works for arbitrarily ordered sums, where one instead considers finite 
Bn Z' Z^^'^ and similarly gets that the sequence of the 

(2.141) Y K.(a:)e''"-*(*«'«'') 

is Cauchy in the desired space because the sequence of the X]aeB„ K«(a;)e*"'^ is Cauchy. □ 
2.6 Error analysis 

The following proposition is used to make precise the notion that <Tp^i{x, 0p) and ^{x, 0p) are approximate 
solutions of (PTM)) and respectively. 

Proposition 2.40. Let R(x, 0) G Hj,^{x, 0) have the property that it is polarized by the elliptic projector Eg 
defined in (j2.23p . i.e. that Ee(i?) — R, and suppose 

(2.142) i?U,=o = 0. 
Then 

(2.143) lim \R{x, 0o + OJ^U • . ■ , + ^Add)k,=^ \E^^-\x.eo) = 
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Proof. It suffices to consider the case R{x, dp) £ H^{x, 9p), Ep(i?) = R, for some p E V U Af, and show 

(2.144) lim |i?(x,0o +^^j,erf)lc,=^Uj-i(,,eo) = 0. 
First, we show 

(2.145) Imi sup \R{x,0o + uip^d)\^^^^\H^-i(^,,^g^-^ = 0. 
Note R has an expansion of the form 

(2.146) R{x,ep)^ aj{x)e''^''. 

Thus, noting R{x,9n) G H^{x,9o) C C{xd,H^~^{x' ,9o)), fixing Xd, we get the norm 

(2.147) |i?(a;,0,)|^._.(^,,^^= J2 E + 

ieZp\o|/3|<s-i 

We let = T- and find 

(2.148) i?(a;,0o+^pCd) = (e"^'"^'-^^^'e*'^''(^'^J')«''aj(x))e*^'^« . 
Then it fohows 

jeZp\o |/3|<s-i 

(2.149) < e-^|i-^.l«^ 5: 5: |9f-a,(-)li.(.0(l + b1)'^^^'-'''^ 

jeZp\o|/3|<s-i 

(2.150) = e-|i-^.l«''|i?(.,0o)|^.-.(.,,„) 
Now we show that, as e tends to zero, 

(2.151) sup \R{x,9o+ui—)\jjs-i,,g.-^0, 
and 

(2.152) sup |i?(a;,0o+^pT^)|^^-i(,,,eo)^O- 

Set /i(a;(j) = R{x' ,Xd,9o) £ H^^{x' ,9o). To see (|2.151l) . observe that the term is bounded by 

(2.153) sup e-^\^-^X-f\h{xd)\l .^^,^,^^< sup ^^(.'..o)' 

which converges to as e tends to zero, since h e C{xd,H'jr^{x' ,9o)) with /i(0) = 0. That (^TS^ holds 
follows from the fact that this term is bounded by 

(2.154) sup e-|i-^.l^|Mx.)|^.-.(,,,,„) < sup |Ma..)lV^(.',.o)^"'""'^'''*' 
which also converges to as e tends to zero. Now we check that 

(2-155) lim \R{x, 9q + uipS^d)\^^=^ \L^xa,H-^{x' A,)) = 0- 

It is not hard to show that 

(2.156) \R{^,e^+^J^)\hi^Xd,H^^ix',e,)) < |e-'""-''''^lLH^(a;,^?o)|i.(.„,fj(.-,eo)) < D^e, 

which completes the proof. □ 
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Now we are ready to show that, as e — >■ 0, the approximate solution u° converges to the exact solution 
Uf: in L°°. This is a corollary of the following theorem. 

Theorem 2.41. For Mo = 2{d + 2) + l and s > 1 + [Mo + ^] let G{x', do) & H^^^ have compact support 
in x' and vanish in t < 0. Let Ue{x,0o) G E^^ be the exact solution to the singular system for < e < co 
given by Theorem 7.1 of /UF, let V° G H^+^ be the profile given by (|1.30[) . and letU^ E £^ be defined by 

(2.157) U\x,eo,^d) ^V'^ix^Oo + u^U. ■ ■ ■ .Oo + ui^i^d). 

Here < To < T is the minimum of the existence times for the quasilinear problems p.36p and (|2.4ip . 
Define 

(2.158) U^{x, 0o) := Z^°(a;, 0o, y ). 
The family lA^ is uniformly bounded in E^^ for < e < eg and satisfies 

(2.159) |J7, - ^ as e 0. 

Upon evaluating our placeholder = ^ in the argument of W G we get the function in E^, 
the same space containing the solutions of the singular system. Indeed, the E^ norm is that of the estimate 
of Proposition 12.431 key to our proof of Theorem I2.41[ and the norm in which we show our approximate 
solution W° is close to the exact solution of the singular system. 

Lemma 2.42. For m > suppose V{x,e) e HJ'+\ EV = V, and setU{x,9o,U) = V{x,9o + uiiU, ■ ■ ■ ,^0 + 
uJ_M^d)- Then, setting U^{x,9o) =U{x,9o, ^), 

(2.160) mE^<\U\£^. 

Proof. The proof is almost identical to the argument used in ^ to prove Lemma 2.25 (b) with Lemma 
2.7. □ 

Proof of Theorem \2.41\ It suffices to prove boundedness of the family in E^^ along with the following 
three statements: 

(a) lim [/" — C/g in E^T^^ uniformly with respect to e G (0, eq] 

n— >-oo ^ 

f2 161) (^) lim Z^,"'" = in E^f^ uniformly with respect to e G (0, eo] 

(c) For each n lim \U^ - p.-i = 0. 

The uniform boundedness of Ue and the [/" in E^^ and that ()2.16ip (a) holds are proved in [11] , Theorem 
7.1, with the use of the iteration scheme (|1.37p and the following linear estimate: 

Proposition 2.43 ([TT], Cor. 7.2). Let s > [Mq + ^] and consider the problem ([07)) . where G G H^^^ 
has compact support and vanishes int<0, and where the right side of (jl.37l) f a) is replaced by € with 
supp T C {t > 0,0 < Xd < E}. Suppose t/" G E^ has compact x— support and that for some K > 0, ei > 
we have 

(2.162) \U:\e^ + \ed,,U:\L^ < K for e e (0, ei]. 

Then there exist constants To{K) and eo{K) < ei such that for < e < eo and T < To we have 

(2.163) \u:'+'\e^+Vt{u:+')s+i,t < CiK,E)VT{\T\E^ + {G)s+i,t) ■ 

The desired boundedness of and the in E^^ follows from the boundedness of V° and the uniform 
boundedness of the V°'" in H^^^, considered with Proposition 12.391 and Lemma [2.421 Similarly, from the 
fact that 

(2.164) lim V"'" = V" in Hf, , 
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we get (|2.16ip fbV 

Now we proceed by proving (|2.161l) (c). Fix S > 0. Noting EV°'" V°'" and EV°'"+i = yo.n+i^ choose 
finite partial trigonometric polynomial sums V^'" and Vp'""*"^ of (resp.) V°'" and V"^""*"^ such that 

(a) EV°'" = V°'" and EV°'"+i = V^"'+\ 

and note that, as a result of (|2.165p fb). we also have 

(2.166) |a,,V°'"+i - a,,Vp"'"+ < '5. 
Observe then, as a consequence of Proposition 12.391 we get 

(2.167) \U°^"-U°'"\£. < C6, < C5, and - d,Mp''^\'-^ < C5, 

where and are obtained from (resp.) V^^" and V°'"+i via the substitution 9 = 6'(6'o,Cd)- The 

induction assumption is 

(2.168) lim|C/:-ZYP'"|^.-. =0. 
With the boundedness of the in , it follows 

(2.169) lim \F{eU^){U:) - F(0)Z^°-"|^.-. = 0. 

Observe that since V°'" and V^^" are invariant under E, Lemma 12.421 applies. Thus, using the estimate 
(|2.160l) . from (|2.169p and (|2.167l) we obtain 

(2.170) \F{eU-){U-) - i^(0)W°;rU.-i <C5 + c{e), 

Tq 

where c(e) — ?> as e — )• 0. 
Now we define 

(2.171) Gp = Z(9,)V°'"+i + A^(Vp"'")9eV°'"+i. 
Then 

- FiO)V^-") =E(L(a,)(V°'"+i - V°^"+i)) + 

(2.172) E(7W(V°'")aeV^^"+i - X(V"'")(99V°^"+^) + 

E(F(0)(V"^"- V;]^"))+i?"+\ 

Thus using (|2.165p fb) and continuity of both E : H"''^ H"'^ and multiphcation from H"'^ x H"'^ i/'*'^, 
we get 

(2.173) lEQp - E{F{0)V"p") - R"+'\h^^ - 0(<5). 
We define the operator 

(2.174) Lo = L{d^) + -i{dcpo)deg + M' {Uli-^)de,. 
We claim 

(2.175) |Lof/r+i - ^^(6t/r)(7r|^.-i < + c(e), 

^0 
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where c(e) — > as e — > 0. This foUows from (|1.37p fa) and 



'To 



0(e) 



(2.176) 



1 7 



1 -7 



-MeU-)f3,de„U:+' - -A,(0)/3,9,„C/r+i + 9„A,(0)[/r/3,9,„C/r+i 



0{e) 



duMo)iu:-u^>;:)Pjde„u, 



rn+l 



^<C\U:-U^>■-\^.-^<cie) + OiS) 

Tq 



"To 



Setting g'p = L{d^)U°-"+' +M'{U^^")de„U^-''+\ since £'(9^0, %)Z^p°'"+' = 0, we get 



(2.177) 



0,n+l _ ^1 



From now on we use \{eo,id) denote evaluation aX — 9{9o,£^d), and |(eo.5j).e to indicate \{eo.id) foUowed 
by evaluation at = It is easy to check Q'^ = Gp\{eo,id)' particularly because V^'", Vp'""*"^ are finite 
polynomials. Thus 



0,n 



0,n 



(2.178) 

(2.179) 

We claim: 
(i) 

(2.180) 
(ii) 

(2.181) 



[ep-F(0)V°^"]|(,„,5,),, 

m{Qv - i^(0)Vp"'") - i?"+i) + + (/ - E)(5p - F{Q)Vl-n]Wu) 



To 



\un\R-+\eo,,dU 
and (iii) there exists Vp G H^^ such that for ~ '^p\o{eo,id)^ 

(2.182) /:'(a,„,a^jz^pi = -[(/ - E){gp - F(o)v°'")]|e(e„,^,). 

(i) follows from ()2.173p followed by application of Proposition 12.391 and then Lemma 12.421 
To see (ii), observe it is clear from our definition of that the continuity of E implies R"^^{x,6) E 

H^^{x, e), and recall from Remark [233] that we have Eei?"+i = and (|2.12ip . By requiring and 

thus (j2.92p . to hold at x^i = 0, we have obtained i?"+^|2,^=o = 0. So Proposition [OOl grants us (I2.18ip . 
Regarding (iii), note that, by Remark 12.81 (|2.182p holds if and only if 

(2.183) £'{deo,di,)U'p=-[{I-E'){gp-F{0)V'p-)]\o^Oo,U, 
for which a sufficient condition is that, where Up = ^p|e(eo,Cd)' 

(2.184) £(9e)Vpi = -(/ - E'){gp - F(0)V°^"). 

Proposition [2l9] guarantees the existence of such Vp, so proof of (iii) is complete. 
Noting 

(2.185) Lo(eWi J = (£'(a«„, 5j jZ^p^), + (L(9.)eWpi). + A^'^D^.o (eZ^p^J, 
it follows from (|2.179p - (|2.182p that 

(2.186) + eUl,) ~ F{0)U^;^\^.^-i < CS + c(e) + K{6)e, 
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where c(e) has been altered, but still tends to zero as e — > 0. It follows from equations (|2.170p . (|2.175l) . and 
(|2.186P that 



|Lo(t/r+i - + dAl,))\E^^-. <C5 + c{e) + K{5)e. 



(2.187) 

Now we claim that we have the following estimates: 



(2.188) 



(a) 



< C(5 + c(e) + ii'((5)e 



(h) |i?(eZ^O-") (C/r+^ - + 6Z^p\J) 1^. < C<5 + c{e) + X(<5)e. 



That (|2.188p fa) holds follows from (|2.187p with the use of estimates similar to (I2.176p . and (j2.188p (b) follows 
from the boundary conditions (|1.36p fb) and (I2.122p (c) together with Proposition l2 .391 and Lemma r2.42l After 
applying the estimate of Proposition 12.431 we get that 



(2.189) 

from which we conclude 
(2.190) 



|C/:+i - (W°;r+i + < C<5 + c(e) + K{5)e. 

To 



< C5 + c{e)+K{5)<i 



This finishes the induction step, completing the proof of the theorem. 



□ 
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